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CHAPTER  1 

INTRODUCTION 


1.1  Background 

Cavitation  nuclei  is  the  general  name  used  to  refer  to  the 
impurities  which  cause  weak  spots  in  a  liquid,  and  thus  prevent  the 
liquid  from  supporting  the  amount  of  tension  which  theory  would  predict. 
It  is  this  inability  of  a  liquid  to  support  a  tension  that  leads  to 
the  phenomenon  of  cavitation. 

Various  forms  of  cavitation  nuclei  have  been  theorized,  but  all 
can  be  considered  in  one  of  two  classes,  depending  on  where  they 
originate.  Stream  nuclei  exist  in  the  free-stream  liquid  in  such  forms 
as  solid  particulates  or  stabilized  microbubbles.  Surface  nuclei 
originate  in  the  surface  of  the  fluid  boundary,  by  means  of  cracks 
and  crevices  in  the  boundary.  A  significant  effort  has  been  devoted 
to  studying  the  nature  of  these  cavitation  nuclei,  and  an  excellent 
review  of  this  research  is  given  by  Holl  in  Reference  [1],  Although 
a  qualitative  relationship  has  been  established  between  the  presence 
of  cavitation  nuclei  and  the  appearance  of  cavitation  during  experi¬ 
ments,  progress  in  making  quantitative  predictions  of  such  relation¬ 
ships  has  been  dependent  upon  the  ability  to  obtain  measurements  of 
cavitation  nuclei  properties  such  as  shape,  size  distribution,  and 
concentration.  Therefore,  much  research  has  also  been  concerned  with 
developing  systems  to  obtain  such  measurements.  Morgan  [2]  gives  a 
summary  of  various  methods  employed  to  measure  cavitation  nuclei. 


In  investigating  the  cavitation  phenomenon  at  ARL/PSU,  progress 
was  hampered  by  the  need  to  measure  the  stream-type  cavitation  nuclei 
distribution  during  cavitation  experiments.  To  remedy  this,  a  system 
to  obtain  such  measurements  was  developed,  patterned  after  the  method 
introduced  by  Keller  [3].  This  method  was  based  on  measuring  the 
scattering  of  a  laser  light  beam  as  nuclei  moved  through  it,  and  was 
thus  called  a  laser  light -scattering  system,  often  abbreviated  as  LSS. 

The  Keller  system  was  chosen  because  it  was  believed  to  have  the 
greatest  potential  for  meeting  the  requirements  deemed  necessary  in  a 
system  for  measuring  cavitation  nuclei  during  experiments.  These 
requirements  arose  because  of  the  peculiar  experimental  conditions 
involved  in  cavitation  research,  and  demand  the  following  capabilities: 

1.  Obtain  measurements  in  the  flow  environment  itself,  yet 
without  disturbing  it:  It  is  a  well-known  fact  that  the 
occurrence  of  cavitation  is  highly  sensitive  to  the  geometry 
of  the  fluid  boundary  and  the  nature  of  the  fluid  flow 
Itself.  For  this  reason,  the  measurement  system  must  not 
affect  the  fluid  flow,  yet  it  must  be  able  to  obtain 
measurements  at  the  very  specific  locations  on  the  fluid 
boundary  where  the  cavitation  occurs. 

2.  Obtain  measurements  on-line:  It  has  also  been  shown  that 
the  cavitation  nuclei  distribution  cannot  be  assumed  to 
remain  constant  with  respect  to  time,  a  phenomenon  which 
has  been  linked  to  the  lack  of  repeatability  in  cavitation 
experiments  [5],  Therefore,  in  order  to  conduct  meaningful 
investigations  into  this  kind  of  phenomenon,  the  measurement 


system  must  give  information  about  the  nuclei  distributions 
which  exist  during  cavitation  experiments,  i.e.,  on-line. 

3.  Separate  measurements  of  solid  particulate  nuclei  from 
gas  bubble  nuclei  in  the  same  flow:  The  necessity  of 

this  requirement  became  apparent  after  a  series  of  experiments 
conducted  jointly  by  the  Applied  Research  Laboratory 
at  Penn  State  University  and  the  California  Institute  of 
Technology.  One  of  the  conclusions  made  from  this  effort 
was  that  differences  in  cavitation  behavior  in  test  facilities 
were  attributable  to  the  variation  in  the  gas  bubble  fraction 
of  the  total  nuclei  population  [6] .  Thus,  it  was  considered 
essential  that  a  nuclei  measurement  system  be  able  to  screen 
out  the  solid  particulate  nuclei  and  examine  only  the  gas 
bubble  distribution. 

4.  Measure  nuclei  in  the  size  range  of  5  to  80  microns  in 
diameter:  Experience  and  theory  have  shown  that  stream- 
type  nuclei  which  affect  cavitation  behavior  are  in  this 
range. 

As  the  Keller  system  was  used  in  some  recent  experiments  [4], 
it  became  clear  that  there  were  still  some  fundamental  inaccuracies  and 
uncertainties  in  the  system’s  performance  that  would  have  to  be  defined 
and  eliminated  before  it  could  be  relied  upon  to  give  the  kind  of  data 
needed  to  conduct  further  meaningful  research  into  the  cavitation 
phenomenon.  For  example,  an  attempt  was  made  by  Keller  to  distinguish 
between  gas  bubbles  and  solid  particulates  by  subtracting  distributions 
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obtained  at  high  pressure  from  those  at  low  pressure,  essentially 
assuming  that  the  gas  bubble  distribution  will  be  eliminated  at  high 
pressure,  leaving  only  the  solid  particulate  distribution.  This  requires 
some  basic  knowledge  of  the  statistics  of  the  flow  which  are  not 
generally  known.  In  addition,  the  pressure  cannot  be  increased  high 
enough  in  most  flow  situations  to  apply  this  method.  In  any  case,  the 
use  of  this  method  reveals  that  the  Keller  system  cannot  distinguish 
bubbles  and  solid  particulates  independently. 

Another  major  inaccuracy  observed  in  the  Keller  method,  as  noted 
in  Reference  [4],  is  the  tendency  for  gas  bubble  concentrations  of  a 
given  radius  to  register  counts  not  only  in  the  distribution  channel 
corresponding  to  that  radius,  but  also  in  all  channels  lower  than 
that  radius.  The  overall  effect  is  that  the  size  distribution  reported 
by  the  system  is  significantly  different  from  the  actual  distribution. 

1.2  Objectives  of  the  Investigation 

It  was  the  intent  of  this  investigation  to  set  forth  the  principles 
involved  in  the  LSS  inaccuracies,  to  describe  methods  whereby  these 
inaccuracies  could  be  minimized  or  eliminated,  and  finally  to  implement 
these  methods  as  modifications  to  the  present  LSS.  This  was  to  achieve 
the  goal  of  providing  a  more  reliable  system  for  the  measurement  of 
cavitation  nuclei. 

In  order  to  accomplish  these  objectives,  the  scope  of  this 
investigation  covered  the  following  issues: 

1.  The  establishment  of  a  criterion  for  determining  the  optimal 
size  of  the  nuclei  measurement  region  in  the  medium. 
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2.  The  development  of  a  system  which  would  be  able 

to  separate  measurements  of  microbubble  nuclei  from 
measurements  of  solid  particulate  nuclei. 

3.  The  utilization  of  a  mathematical  model  to  correct  measured 
size  distributions. 

As  an  additional  issue  in  the  scope  of  this  investigation,  a  method  for 
correcting  size  distributions  by  means  of  direct  measurements  of  the 
intensity  distribution  in  the  measurement  region  was  developed. 

The  procedure  of  this  investigation  was  to  obtain  the  needed 
criteria  and  mathematical  models  primarily  from  existing  theories  and 
known  principles  of  the  system's  operation,  and  then  to  conduct 
experiments  that  would  verify  those  theories  .  The  plan  involved 
taking  a  large  number  of  measurements  of  many  different  nuclei 
distributions,  these  using  the  original  Keller  system,  a  modified  LSS 
system,  and  a  holography  system  which  would  provide  a  calibration  of 
the  true  nuclei  distributions.  Data  from  these  experiments  were  then 
compared  with  predicted  results  in  determining  the  optimal  measurement 
region  (or  Probe  Volume)  size,  and  similarly  for  validating  the  method 
of  separating  bubble  nuclei  and  solid  particulate  nuclei.  This  same 
data  base  was  also  used  to  both  implement  and  validate  the  mathematical 
model  to  correct  measured  bubble  size  distributions. 

The  latter  part  of  the  experiments  was  to  obtain  measurements  of 
the  probe  volume  intensity  distributions  for  several  conditions,  and 
then  to  correlate  these  measurements  with  measurements  of  the  bubble 
size  distribution  taken  simultaneously.  This  data  base  provided  the 
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means  for  implementing  and  verifying  the  method  of  correcting  bubble 
size  distributions  through  direct  measurement  of  Probe  Volume  intensity 
distributions . 


i 


\ 


i 


i 

1 


7 


CHAPTER  2 

DESCRIPTION  OF  THE  LIGHT-SCATTERING  SYSTEM 

2.1  The  Keller  System 

2.1.1  Principle  of  Operation.  The  governing  principle  of  this 
measurement  system  is  dependent  on  the  size  regime  of  the  nuclei 
being  measured.  Measurements  in  the  micron  size  range  dictate  that 
the  Mie  scattering  theory  be  used,  since  this  is  the  theory  relevant 
for  nuclei  whose  sizes  are  on  the  same  order  of  magnitude  as  the  wave 
lengths  of  the  laser  light  which  illuminates  them.  However,  for  the 
measurement  of  the  larger  microbubbles,  geometric  theories  could  be 
utilized. 

Mie  theory  predicts  that  a  monochromatic,  coherent,  plane-wave 
light  beam  will  be  scattered  in  all  directions  by  an  arbitrary  particle 
and  tha  the  intensity  of  that  scattered  light  will  depend  strongly 
on  the  size,  shape,  and  refractive  index  of  the  particle.  The  intensity 
will  also  vary  greatly  with  the  solid  angle  into  which  the  light  is 
scattered.  For  spherical  particles,  Mie  theory  predicts  that  in  certain 
directions  perpendicular  to  the  illuminating  beam  direction,  the 
scattered  light  intensity  will  increase  monotonically  for  a  limited 
particle  size  range  and  a  given  refractive  index.  It  is  this  characteristic 
of  the  scattering  theory  which  forms  the  basis  of  the  measurement 
system,  as  shown  in  Figure  1.  A  derivation  of  this  relationship  is 
given  in  Appendix  A. 


z 


DIRECTION 

Figure  lb.  Scattering  Coordinate  System  for  the  ARL/PSU  LSS. 


Figure  2.  Simplified  Representation  of  Probe  Volume  (PV).  Dimensions. 
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2.1.2  Method  of  Operation.  The  actual  cavitation  nuclei  measure¬ 
ment  system,  which  is  based  on  the  Keller  method,  can  be  thought  of 
as  being  composed  of  three  primary  subsystems: 

1.  Sensing  subsystem  -  The  purpose  of  the  sensing  subsystem 
is  to  construct  a  monochromatic,  coherent,  plane-wave 
light  beam  of  known  dimensions  and  intensity  at  a  location 
in  the  flow  where  the  nuclei  distribution  is  desired  to  be 
measured  (this  location  will  hereafter  be  referred  to  as 
the  Probe  Volume  (PV)).  This  purpose  is  achieved  in  the 
Keller  system  by  passing  the  beam  from  a  laser  light  source 
through  the  optical  system,  such  as  is  shown  in  Figure  3. 

The  primary  issue  in  the  sensing  subsystem  is  to  obtain  a 
beam  at  the  PV  that  has  optimal  intensity  and  physical 
dimensions  so  as  to  minimize  statistical  errors  encountered 
when  sampling  large  numbers  of  nuclei.  The  nature  of 
these  errors  will  be  discussed  further  in  Chapter  3. 

2.  Measuring  subsystem  -  The  purpose  of  the  measuring  subsystem 
is  to  record  the  amount  of  laser  light  which  is  scattered 

by  a  nucleus  as  it  passes  through  the  laser  beam 

at  the  PV  (the  laser  beam  at  the  PV  which  illuminates  the 

nuclei  will  often  be  referred  to  as  the  incident  beam) . 

The  Keller  system  accomplishes  this  by  using  a  photo¬ 
sensitive  detector,  or  photomultiplier  tube  (abbreviated 
as  PM-tube),  to  detect  the  intensity  of  scattered  light. 

As  seen  in  Figure  3,  the  scattered  light  is  collected  by 
a  set  of  lenses  whose  diameter  and  focal  length  determine 
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Figure  3.  Schematic  of  the  Keller  LSS . 
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a  fixed  solid  angle  of  scattered  light  to  be  collected. 
This  light  is  then  focused  onto  a  slit  aperture,  which 
restricts  the  length  of  the  laser  beam  monitored  by  the 
PM-tube.  The  PM-tube  then  converts  the  light  scattered 
by  individual  nuclei  into  equivalent  electrical  signals. 
The  slit  aperture  restricts  the  amount  of  light  entering 
the  PM-tube  in  order  to  set  the  third  dimension  of  the 
laser  light  PV,  the  other  two  dimensions  are  set  by  the 
cross-sectional  area  of  the  incident  beam  (see  Figure  2) . 
The  primary  issues  in  this  subsystem  are  the  providing  of 
a  guarantee  that  the  PM-tube  is  at  the  beam-perpendicular 
location  shown  to  be  optimal  by  the  Mie  theory,  and 
optimization  of  the  size  of  the  slit  aperture  so  as  to 
minimize  the  statistical  errors  mentioned  in  the  sensing 
subsystem. 

3.  Analyzing  subsystem  -  The  purpose  of  th<L«  vjbsystma  is 

to  accurately  interpret  the  recordings  of  scattered  laser 
light  obtained  over  a  fixed  time  interval,  transforming 
the  electrical  signals  from  the  measuring  system  into  a 
histogram  representing  the  nuclei  size  distribution  which 
was  present  during  that  time  interval.  This  is  achieved 
in  the  Keller  system  by  passing  the  electrical  signals 
from  the  PM-tube  into  a  device  called  a  Pulse-Height 
Analyzer  (PHA)  to  obtain  the  distribution  histogram. 

The  PHA  is  a  device  which  can  be  thought  of  as  the 
electronic  equivalent  of  the  scattering  relationship 
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shown  in  Figure  1,  with  the  difference  being  that  the 
relationship  between  nuclei  size  and  scattered  light 
intensity  is  not  continuous  but  is  coarse-grained  into 
discrete  intervals,  called  PHA  channels.  The  number  and 
relative  width  of  the  PHA  channels  are  arbitrary,  being 
governed  primarily  by  the  degree  of  measurement  accuracy 
desired.  The  primary  issue  in  the  analyzing  subsystem 
is  to  interpret  a  large  enough  data  base  so  as  to  minimize 
statistical  errors,  yet  do  it  in  a  small  enough  time 
interval  to  still  qualify  the  light-scattering  system  as 
an  on-line  measurement  system. 

2.2  Identifying  System  Inaccuracies 

In  this  investigation  the  inaccuracies  under  consideration  have 
their  origin  in  each  of  the  three  subsystems  detailed  in  Section  2.1.2. 

In  the  sensing  subsystem,  a  need  for  a  clearer  definition  of  the  size 
and  shape  of  the  PV  was  observed,  so  as  to  better  understand  the 
amount  of  error  inherent  in  the  nuclei  distribution  measurements. 

The  confidence  in  the  system’s  measurements  is  in  direct  proportion 
to  the  precision  of  the  system's  error  estimate;  thus,  such  a  definition 
of  the  PV  and  the  errors  associated  with  it  was  deemed  essential  for 
a  reliable  measurement  system. 

The  inability  of  the  system  to  meet  the  requirement  of  distinguishing 
between  microbubbles  and  solid  particulates  was  traced  back  to  an  under¬ 
utilization  of  the  scattered  light  signals  as  recorded  by  the  measuring 
system. 
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A  major  inaccuracy  was  observed  in  the  analyzing  subsystem, 
whereby  the  final  nuclei  distribution  reported  by  the  PHA  was  consistently 
found  to  be  significantly  different  from  the  actual  distribution. 
Experiments  done  in  Reference  [7]  showed  that  individual  scattering 
from  nuclei  does  repond  in  the  manner  predicted  by  the  Mie  theory, 
so  it  was  concluded  that  the  problem  must  be  limited  to  the  incident 
beam  character  itself  in  the  sensing  subsystem. 


i 
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CHAPTER  3 

PRINCIPLES  OF  INACCURACIES  AND  PROPOSED  SOLUTIONS 

3 .1  The  Criteria  for  the  Probe  Volume 

3.1.1  Nature  of  Potential  Errors.  As  was  discussed,  the  size 
of  the  laser  light  Probe  Volume  (PV)  is  a  source  of  potential  inaccuracy 
in  the  light-scattering  system  (LSS)  .  This  inaccuracy  arises  because 
not  all  microbubbles  of  the  same  size  which  enter  the  PV  will  scatter 
the  same  amount  of  light.  There  are  two  mechanisms  which  can  cause 
this  departure  from  predicted  scattering  behavior:  coincidence  errors 
and  illumination  errors. 

Measurement  inaccuracy  due  to  coincidence  errors  is  usually 
caused  by  the  smaller  nuclei  in  the  range  of  sizes  being  measured. 

If  the  PV  is  large  enough,  it  is  possible  for  more  than  one  scatterer 
to  be  illuminated  simultaneously  by  the  laser,  yet  only  a  single 
composite  light-intensity  measurement  is  recorded  by  the  photo¬ 
multiplier  tube  (PM-tube) .  The  resulting  measurement  goes  on  to  be 
interpreted  as  a  single  large  nucleus,  instead  of  several  small  ones, 
and  an  incorrect  picture  of  the  size  distribution  is  obtained. 

Inaccuracy  due  to  illumination  errors  is  caused  by  the  scatterer 
in  the  PV  not  being  illuminated  by  the  peak  intensity  of  the  laser 
light.  If  the  scatterer  is  large  enough,  it  is  also  possible  for  it 
to  pass  through  an  "edge"  of  the  PV  so  that  only  part  is  illuminated. 
Both  of  these  situations  cause  the  LSS  to  interpret  the  corresponding 
scatterer  as  being  smaller  than  it  actually  is,  again  leading  to  an 
incorrect  picture  of  the  size  distribution. 
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3.1.2  Coincidence  Error  Criterion.  The  objective  of  a  coincidence 
error  criterion  is  to  establish  a  formula  whereby  the  extent  of 
inaccuracy  in  measurements  due  to  coincidence  effects  can  be  calculated 
for  a  given  set  of  experimental  conditions.  Any  realistic  probability 
model  of  the  interaction  of  a  scatterer  with  the  PV  must  be  based  on  the 
Poisson  probability  distribution,  which  gives  the  probability  of  k 
scatterers  being  simultaneously  in  the  PV  as 

k  -x 

P(k,x)  =  x— 7—  (3.1) 

where  x  =  C  T 

C  =  average  total  rate  of  scatterers  passing  through  the  PV 
X  =  PM-tube  resolving  time  (i.e.,  average  length  of  time  to 
pass  through  the  PV) . 

It  can  be  seen  from  Equation  (3.1)  that  the  probability  of  one  or  more 
scatterers  being  simultaneously  in  the  PV  (in  other  words,  the 
probability  of  interference)  is  given  by 

Pint  -  1  -  P(o,x).  (3.2) 

Analyses  to  estimate  P^nt  have  been  done  with  varying  degrees  of 
depth  by  several  authors  [7,8,9,10,11],  In  most  cases,  P  was  set 
at  some  specified  value  acceptable  to  the  experimenter  and  was  then 
related  to  other  experimental  parameters.  The  final  formula  involved 
a  choice  between  specifying  either  the  size  of  the  control  volume,  or 
the  maximum  total  concentration  to  be  measured. 

The  most  rigorous  treatment  was  considered  to  be  that  done  by 
Holve  and  Self  [10],  Their  analysis  began  with  the  Poisson  probability 
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distribution  in  the  manner  of  Equations  (3.1)  and  (3.2),  and  then 
proceeded  in  terms  of  important  experimental  quantities  such  as  the 
size  range  able  to  be  measured  and  the  total  nuclei  concentration,  as 
well  as  an  estimate  of  the  overall  size  concentration  distribution. 
The  resulting  formula  for  estimating  coincidence  errors  was 


nTVm 

3  =  1JL. 

int  (q-1) 


lo8(d«*/d»in) 


(3.3) 


where  P .  =  acceptable  probability  of  coincidence 

V  =  probe  volume  size 
m 

$  m  ratio  of  largest  to  smallest  size  in  the  range 

max  min 

that  the  PHA  is  set  up  to  measure 
nT  =  total  nuclei  concentration 
q  *  size  concentration  distribution  parameter. 

The  variable  q  in  Equation  (3.3)  comes  from  the  assumption  that  the 
form  of  the  size  concentration  distribution  can  be  estimated  as 


N(d)  -  Ad~q,  (3.4) 

where  d  is  the  scatterer  size,  A  is  a  constant,  and  q  is  always  positive. 

Due  to  its  rigorous  derivation,  Equation  (3.3)  can  be  viewed  as  the 
general  form  of  the  error  estimation  formula,  of  which  the  estimation 
criterion  of  other  investigators  are  special  cases.  For  instance,  in 
Reference  [11]  Keller’s  formula  for  estimating  coincidence  errors  is 


int 


(3.5) 
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In  both  References  [7]  and  [8],  the  conclusion  is  that  the  PV 
size  which  minimizes  coincidence  errors  is  given  by 


1 


n_V 
T  m 


(3.6) 


In  the  present  investigation,  Equation  (3.3)  was  used  to  estimate 
coincidence  errors,  and  the  experimental  conditions  involved  resulted 
in  the  following  formula: 


P 


int 


n^V 
T  m 

1.43 


(3.7) 


The  actual  conditions  which  gave  rise  to  Equation  (3.7) ,  together  with 
an  analysis  of  the  resulting  experimental  data,  are  given  in  Section  5.1. 


3.1.3  Illumination  Error  Criterion.  In  treating  the  inaccuracies 
associated  with  improper  scatterer  illumination,  one  must  deal  with 
the  problem  of  non-uniformity  in  the  laser  beam  intensity  distribution 
Itself.  In  Appendix  A  of  Reference  [11],  Keller  attempts  such  a  treatment 
by  postulating  regions  of  uncertainty  surrounding  the  control  volume. 

The  assumption  is  made  that  the  ratio  of  the  Incorrect  measurements 
(due  to  illumination  errors)  to  the  total  number  of  measurements  is 
equal  to  the  ratio  of  this  uncertainty  area  to  the  total  area  of  the  PV, 
normal  to  the  flow.  This  model  of  the  PV  and  the  resulting  error 
estimation  formula  are  given  in  Figure  4. 

However,  this  model  assumes  a  uniform  laser  light  intensity  in 
the  region  of  the  PV  which  is  not  part  of  the  uncertainty  region,  an 
assumption  which  Keller  himself  shows  does  not  represent  physical 
reality  [11].  In  fact,  the  laser  beam  has  a  Gaussian  Intensity 
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distribution,  and  this  occurs  only  under  ideal  experimental  conditions. 

It  was  also  observed  in  this  investigation  that  the  Probe  Volume 
illumination  intensity  could  not  be  assumed  to  be  uniform,  or  even 
Gaussian.  Therefore,  it  was  concluded  that  a  fundamentally  different 
analysis  must  be  used  to  estimate  the  effect  of  illumination  errors, 
one  which  would  take  into  account  the  non-uniformity  of  the  laser  beam. 
Such  an  analysis  is  detailed  in  Section  3.3. 

3.2  The  Modified  Measuring  Subsystem 

3.2.1  Principle  of  the  Dual-Detector  Operation.  As  mentioned 
previously  in  the  system  description,  the  inability  of  the  LSS  to 
distinguish  between  microbubble  nuclei  and  solid  particulate  nuclei 
is  due  to  an  under-utilization  of  the  full  Mie-scattering  properties 
of  the  nuclei.  The  potential  for  this  distinguishability  comes  through 
the  solution  to  the  transverse  angular  solution  of  the  full  Mie  scattering 
equations  (see  Equation  (A. 4)) .  This  solution  has  the  form 

0(0)  -  P^m)(cos0),  (3.8) 

where  are  the  Legendre  polynomials.  The  fact  that  the  argument 

of  the  Legendre  polynomials  is  cos0  indicates  that  the  solution  is 
symmetric  at  about  0=0.  Therefore,  the  intensity  of  scattered  light 
at  0  *  90°  will  be  the  same  as  at  0  =  -90°,  but  only  if  the  nuclei  are 
ideal  Mie  scatterers,  that  is,  spherically  symmetric. 

This  property  of  the  scattering  equations  allows  one  to  take 
advantage  of  the  difference  in  geometry  between  microbubble  and  solid 
particulates.  The  microbubbles,  being  spherical,  will  cause  symmetric 
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light-scattering,  whereas  the  particles  are  in  general  randomly  shaped 
and,  therefore,  will  scatter  light  randomly.  In  principle,  the  way  to 
separate  measurements  of  bubbles  from  measurements  of  solid  particulate 
is  to  compare  the  scattered  light  intensity  from  both  0  =  90°  and 
0  =  -  90°.  In  practice,  the  issue  is  the  establishing  of  a  criterion 
for  the  cutoff  asymmetry  of  scattering:  That  is,  how  equal  must  the 
scattered  intensities  at  these  two  locations  be  in  order  for  the 
scatterers  to  be  considered  as  microbubbles;  or  similarly,  how  unequal 
must  they  be  in  order  to  be  considered  as  solid  particulates? 

3.2.2  A  Theory  for  Scattering  by  Non-Symmetr ic  Particles.  This 
question  has  been  addressed  in  a  recent  analysis  by  Kohler  and  Billet  [12] 
where  the  asymmetry  ratio  of  a  general  randomly-shaped  particle  is 
calculated  as  a  function  of  its  orientation  (with  respect  to  the  incident 
beam  axis),  and  as  a  function  of  its  axial  ratio  and  refractive  index. 

The  asymmetry  ratio,  R,  is  defined  as  the  ratio  of  the  smaller  magnitude 
scattered  intensity  to  the  larger.  The  axial  ratio  is  defined  as  the 
ratio  of  the  particles  semi-minor  axis  to  its  major  axis,  and  comes 
from  considering  the  general  randomly-shaped  particle  to  be  an  ellipsoid 
or  revolution  or  a  prolate  spheroid. 

In  this  analysis,  scattering  calculations  were  made  for  many 
positions  within  the  solid  angle  subtended  by  the  receiving  optics. 

These  calculations  were  performed  for  a  number  of  different  particle 
orientations  to  describe  the  scattering  in  both  directions  (0  *  +  90°) 
as  functions  of  orientation  angles,  then  the  ratios  of  these  scattered 
intensities  resulted  in  R.  Integrating  R  over  all  possible 
orientations  of  the  random  particle  gave  an  average  intensity  ratio,  R. 
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This  same  integration  procedure  was  used  to  compute  the  standard 
deviation  (SD)  from  the  average  intensity  ratio.  In  analyzing  the 
characteristics  of  R,  they  showed  that  it  is  nearly  independent  of 
both  the  particle  size  and  refractive  index.  The  critical  factor 
was  the  axial  ratio  of  the  particle.  Thus,  the  values  of  R  and  SD 
computed  at  one  axial  ratio  would  be  general  for  all  sizes  and  refractive 
indices. 

Choosing  an  axial  ratio  of  1.76,  Kohler  and  Billet  computed  an 
R  of  0.45  and  a  SD  of  0.33.  Assuming  that  the  random  orientations 
were  distributed  normally,  the  normal  distribution  then  predicted 
that  84%  of  the  orientations  would  have  a  value  of  R  less  than  R  +  SD. 
After  adjusting  the  value  of  SD  to  account  for  computational  errors, 
this  R  value  was  given  as  R  +  SD  =  0.75,  and  was  taken  to  be  the 
cutoff  asymmetry  ratio,  R^. 

This  analysis  thus  provided  the  theoretical  basis  for  separating 
measurements  of  microbubbles  and  measurements  of  solid  particulate 
in  a  laser  light-scattering  system. 

3.2.3  Method  of  Dual-Detector  Operation.  The  schematic  diagram 
for  a  modified  LSS  which  has  this  microbubble  distinguishing  capability 
is  shown  in  Figure  5.  The  difference  from  the  original  system  shown  in 
Figure  3  is  that  now  there  are  two  photo-detecting  devices,  one  at 
0  *  90°  and  the  other  at  9  ■  -90°.  Each  has  its  own  set  of  receiving 
optics  and  its  own  signal  conditioner,  and  the  output  of  these  signal 
conditioners  then  go  into  a  new  electronic  device  called  the 
Pulse-Height  Comparator  (PHC) .  The  PHC  performs  the  comparison  of  the 
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scattered  intensities  (now  in  the  form  of  electronic  signals)  described 
in  the  previous  analysis.  If  the  compared  signals  give  a  value  of  R 
which  is  less  than  R^,  then  the  signals  are  rejected  and  no  further 
signal  processing  occurs.  If  the  compared  signals'  intensity  ratio 
is  greater  than  or  equal  to  R^,  then  the  signals  are  passed  on 
to  the  PHA  to  be  sorted  into  the  appropriate  size  channels. 

In  this  investigation,  experiments  were  conducted  to  test  the 
validity  of  the  cutoff  asymmetry  ratio,  R^,  predicted  by  Kohler  and 
Billet,  and  in  so  doing  prove  the  effectiveness  of  the  dual-detector 
system  as  a  means  of  distinguishing  microbubbles  and  solid  particulate 
cavitation  nuclei.  The  results  of  these  experiments  are  given  in 
Section  5.2. 


3.2.4  Some  Implications  of  This  Method  for  Solving  Control 
Volume  Inaccuracies.  An  important  additional  effect  of  the  dual¬ 
detector  scattered  light  comparison  method  relates  to  the  problem  of 
illumination  errors  described  in  Section  3.1.  It  was  mentioned  that 
one  cause  of  this  kind  of  error  was  large  nuclei  (l.e.,  nuclei  whose 
size  is  of  the  same  order  of  magnitude  as  the  CV)  passing  through 
an  "edge"  of  the  control  volume  so  that  only  part  is  illuminated. 

This  situation  is  a  departure  from  an  ideal  Mie  scatterer  because  in 
order  for  nuclei  to  scatter  light  symmetrically  they  must  be  illuminated 
by  a  uniform  plane  wave.  Therefore,  because  the  nucleus  is  non-unif ormly 
illuminated,  it  will  scatter  light  asymmetrically,  and  thus  be  rejected 
as  a  solid  particle  nucleus  by  the  PHC. 

So,  in  addition  to  providing  a  nuclei-distinguishing  capability  to 
the  LSS,  the  dual-detector  system  also  eliminates  the  need  for  an 


analysis  to  account  for  illumination  errors  caused  by  nuclei  crossing 
an  edge  of  the  PV. 


3.3  The  Modified  Analyzing  Subsystem 

3.3.1  Principle  of  the  Nuclei  Counter  Inaccuracy.  Another  major 
problem  observed  in  the  LSS  was  the  inaccurate  classification  of 
nuclei  according  to  size  and  concentration,  as  reported  by  the  Pulse- 
Height  Analyzer  (PHA) .  This  inaccuracy  is  demonstrated  clearly 
in  Figures  6a  and  6b,  which  show  the  size  distribution  from  a  typical 
nuclei  measurement  experiment  as  given  by  the  methods  of  laser 
light-scattering  and  holography.  The  problem  does  not  lie  with  the  Mie 
theory  calibration  relationship  itself,  nor  does  the  problem  appear 
to  be  caused  by  the  geometry  of  the  optics  used,  since  the  same 
problem  was  documented  by  researchers  using  a  near-forward  scattering 
system  [10].  These  same  researchers  concluded  that  the  problem  was 
caused  by  a  non-uniform  intensity  in  the  volume  of  laser  light  used  to 
measure  the  nuclei  (i.e.,  the  PV) . 

As  mentioned  in  Section  3.1,  such  non-uniformity  in  the  PV 
illumination  intensity  was  clearly  observed  in  this  investigation. 
Figures  12a  and  12b  show  a  survey  of  the  illuminating  laser  beam  (using 
a  30  ym  aperture  attached  to  a  photodetector)  made  before  the  beam 
enters  the  water  tunnel  test  section,  and  then  another  survey  of  the 
beam  inside  the  test  section.  The  obvious  conclusion  to  be  drawn  is 
that  the  intensity  in  the  PV  cannot  be  assumed  to  be  uniform,  even  if 
the  incident  beam  is  made  uniform. 


Distribution  Given  by  the  LSS  for  the  70-pm  Experiment. 
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Figure  6b.  Nuclei  Concentration  Distribution  Given  by  Holography  for  the  70-Utn  Experiment 
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3.3.2  Method  of  Matrix  Inversion.  To  deal  with  the  nor.-uniform 
intensity  problem,  the  researchers  in  Reference  [10]  developed  a 
mathematical  method  to  correct  the  count  rate  data  given  by  their 
particle  counter.  The  method  is  called  an  inversion  scheme,  and  the 
technique  used  to  deal  with  inaccuracies  in  the  nuclei  counter  (i.e., 

PHA)  was  patterned  after  this  method.  The  basis  of  the  technique  is 
to  consider  that  the  governing  Mie  scattering  equation  has  a  variable 
incident  beam  intensity  instead  of  a  constant  one  (see  Equation  (B.l)). 
The  effect  is  that  the  equation  used  to  convert  count  rate  data  into 
concentration  data  is  changed  in  the  following  way.  The  actual  data 
histogram  given  by  the  PHA  is  in  terms  of  count  rate,  that  is,  the 
number  of  nuclei  counted  per  unit  time.  The  physically  relevant 
quantity  is  the  nuclei  concentration.  These  quantities  were  considered 
to  be  related  by  a  constant  as 

C  =  USNi,  (3.9) 

where  =  count  rate  of  nuclei  in  channel  i,  as  reported  by  the  PHA 

U  =  speed  of  nuclei  through  the  control  volume 
S  =  area  of  the  control  volume  normal  to  the  flow  direction 
=  concentration  of  nuclei  of  size  category  i. 

The  effect  of  a  variable  incident  beam  intensity  on  Equation  (3.9) 
is  to  change  it  to  a  matrix  equation  expressed  as 

Ci  -  US^Nj,  (3.10) 

where  S . ,  *  matrix  of  area  elements  which  transforms  the  vector  of 
ij 

concentration  data  into  the  vector  of  count  rate  data. 
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In  order  to  have  a  method  of  transforming  count  rate  data  into 
concentration  data,  which  is  what  is  needed  for  the  nuclei  counter. 
Equation  (3.10)  is  solved  for  N.,  yielding 


i' 


(3.11) 


Physically,  Equation  (3.11)  states  that  the  nuclei  concentration 
of  a  given  size  category  is  not  merely  the  corresponding  count  rate 
from  that  category,  but  is  actually  a  weighted  sum  of  the  count  rates 
from  all  channels  in  the  nuclei  counter,  where  the  weighting  factors 
are  the  inverted  elements  of  the  area  matrix,  S...  The  variation  of 
the  elements  in  this  matrix  reflects  the  variation  of  the  illuminating 
beam  intensity  in  the  PV;  therefore,  the  matrix  should  be  constant  for 
a  given  laser  and  optical  arrangement,  so  that  theoretically  the 
elements  need  only  be  calculated  once  for  a  given  arrangement.  This 
also  implies  that  this  classification  inaccuracy  phenomenon  should  be 
seen  more  clearly  as  the  number  of  nuclei  measured  by  the  system  becomes 
statistically  large,  since  in  time  these  nuclei  will  sample  all  regions 
of  the  PV  and  thus  experience  the  full  extent  of  its  non -uniformities. 

A  more  detailed  summary  of  the  inversion  scheme  developed  by  Holve 

and  Self  [10]  that  formed  the  basis  for  this  analysis  is  given  in  Appendix  B. 

In  this  investigation,  a  series  of  calibration  experiments  were 
carried  out  where  U,  C^,  and  were  measured  under  many  different 
conditions  in  order  to  solve  for  the  elements  of  the  area  matrix, 

Sjj,  so  that  this  matrix  could  be  used  according  to  Equation  (3.11) 
to  correct  the  count  rate  histograms  obtained  with  the  light-scattering 


system. 
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CHAPTER  4 

EXPERIMENTS  TO  TEST  PROPOSED  SOLUTIONS 
4 .1  Experiments  Involving  Nuclei  Distributions 

4.1.1  Experimental  Objectives  and  Methods.  The  objectives  of 
this  series  of  experiments  were  to  determine  the  cutoff  level  for  the 
Pulse-Height  Analyzer  (PHA)  for  the  separation  of  microbubble  and  solid 
particulate  nuclei  and  to  determine  the  elements  of  the  area  matrix, 

S_.  The  definition  of  is  the  key  to  the  operation  of  the 

inversion  scheme  for  correcting  some  inaccuracies. 

The  operating  parameter  of  the  Pulse-Height  Comparator  (PHC)  is 
the  cutoff  rejection  ratio,  Rc>  defined  in  Section  3.2.2.  Defining 
this  parameter  requires  not  only  an  estimate  of  the  scattering  charac¬ 
teristics  of  non-symmetric  particles,  but  also  an  estimate  of  the 
asymmetry  of  sphere  scattering  under  actual  test  conditions.  The  method 
used  to  obtain  those  estimates  was  to  program  a  variable  cutoff 
rejection  ratio  in  the  PHC.  Estimates  of  sphere-scattering  asymmetry 
were  then  obtained  by  measuring  the  number  of  spheres  rejected  as  a 
function  of  for  several  different  sphere  size  monodispersions. 
Similarly,  estimates  of  non-symmetric  particle  scattering  were  obtained 
by  measuring  the  total  number  of  rejected  scatterers  (both  spheres  and 
particles)  for  a  mixture  of  microspheres  and  non-symmetric  particles. 

In  order  to  determine  the  elements  of  the  area  matrix,  , 
measurements  of  count  rate,  C^,  and  concentration,  N^,  were  made  of 
single-size  concentrations  (or  monodispersions)  of  spheres.  The  goal 
was  to  measure  monodispersions  representing  each  of  the  10  size 


categories  in  the  nuclei  classification  system.  Such  a  set  of 


measurements  would  allow  the  elements  of  S..  to  be  solved  for  each 

ij 

of  the  10  columns  in  the  matrix.  The  form  of  this  matrix  would  then 
be  10  by  10,  and  would  be  upper-triangular  (i.e.,  all  elements 
below  the  main  diagonal  are  zero)  .  The  upper-triangular  form  is  a 
result  of  the  fact  that  a  given  scatterer  will  not  appear  larger 
than  its  Mie-predicted  size,  since  this  predicted  size  is  based  on 
the  scatterer  passing  through  the  peak  intensity  portion  of  the  PV. 
Appendix  B  gives  a  more  complete  derivation  of  the  equations  governing 
this  experimental  procedure. 


4.1.2  Description  of  Equipment.  Figure  7  gives  a  simplified 
schematic  diagram  of  the  overall  experimental  set-up  used.  This 
set-up  has  three  independent  subsystems: 

1.  Korad  Holography  System  (KHS)  -  The  holography  system 

was  used  to  calibrate  the  light- scattering  system  (LSS) . 

As  the  LSS  was  measuring  the  concentration  of  a  given 
size  category,  a  hologram  of  a  representative  volume 
of  the  flow  during  the  period  of  the  LSS  measurements  was 
obtained.  Reconstruction  of  a  magnified  image  of  the 
hologram  allowed  the  scatterers  to  be  actually  visualized, 


sized  and  counted  to  give  the  concentration  measurements. 


Nj .  Figure  8  shows  a  schematic  diagram  of  the  equipment 
used  to  reconstruct  the  holograms  and  measure  the 


nuclei  concentrations. 

« 

2.  Water  Tunnel  and  Pressuring  System  -  This  system  provided 
the  means  for  controlling  the  nuclei  population  in  the 


Figure  7.  Simplified  Schematic  of  Experimental  System. 
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Figure  8.  Schematic  of  the  Hologram  Reconstruction  System, 
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experiments  and  is  shown  in  Figure  9.  The  water  tunnel 
itself  was  designed  and  built  especially  for  this 
investigation.  It  was  constructed  of  plexiglas  with  a 
fiberglas  casing  and  had  a  3-inch  by  3-inch  square 
cross-section,  except  for  the  test  section  which  was 
1  inch  by  3  inches  and  6  inches  in  length.  Power  was 
provided  by  a  2  hp  centrifugal  pump,  made  of  cast  iron  and 
treated  with  an  epoxy  sealant  to  inhibit  rust  formation. 
The  pump  motor  was  rated  at  3425  rpm,  but  a  frequency 
inverter  was  connected  to  the  motor  to  obtain  speed 
variation. 

Purity  of  the  water  sample  used  in  experiments  was 
achieved  in  the  following  way.  Distilled  water  (with  a 
filtration  purity  of  2  pm)  was  used  as  the  water  sample 
and  was  contained  in  a  30-gallon  sealed  drum.  A  rust 
inhibitor  (nickel  chromate)  was  added  to  suppress  the 
formation  of  rust  in  the  pump  machinery  and,  therefore, 
minimize  the  contribution  of  rust  particles  to  the  back¬ 
ground  nuclei  concentration.  The  water  in  this  container 
was  then  subjected  to  a  vacuum  of  approximately  5  psia  for 
a  period  of  between  12  and  36  hours  to  remove  as  much  of 
the  air  dissolved  in  the  sample  as  possible.  During  the 
operation  of  the  experiment,  the  water  tunnel  was 
pressurized  to  a  maximum  of  8  psia  in  order  to  suppress 
the  generation  of  bubbles  from  cavitation  in  the  pump 
machinery. 


Schematic  of  the  Water  Tunnel  and  Pressuring  System. 


3.  Light-Scattering  System  (LSS)  -  The  LSS  used  in  the 


experiments  was  different  from  that  detailed  in  Chapter  2 
only  in  that  both  the  old  and  new  PHA's  were  used.  The 
old  PHA,  which  interpreted  data  from  only  one  PM-tube,  was 
used  to  calibrate  the  new,  two-PM-tube  PHA;  that  is, 
the  old  PHA  was  used  to  make  sure  that  the  two  PM-tubes, 
which  feed  into  the  new  PHA,  had  equal  amounts  of  amplifier 
gain.  This  was  important  to  assure  that  both  PM’s  were 
indeed  measuring  the  same  signals  before  comparisons  of 
differing  signals  in  the  PHC  could  be  made  with  confidence. 
The  matching  procedure  was  to  feed  the  output  of  each 
detector  into  an  older  model  nuclei  classification  system, 
which  was  based  on  only  one  output.  The  counting  times 
of  each  detector  would  be  noted  as  they  measured  an 
actual  steady  flow;  then  the  electronic  gain  on  the 
amplifiers  associated  with  the  detectors  would  be  adjusted 
until  both  detectors  would  consistently  give  the  same 
counting  rate.  One  detector  was  the  one  with  which  the 
calibration  relationship  between  nuclei  size  and  detector 
voltage  was  determined,  then  the  second  detector's 
amplifier  gain  was  matched  to  this  one.  Figure  10  shows 
a  schematic  diagram  of  the  equipment  used  in  the 
experimental  version  of  the  LSS. 


4.1.3  Procedure  for  the  Experiments.  The  following  is  the 
chronological  procedure  used  in  the  experiments  at  the  Materials 
Research  Laboratory  of  The  Pennsylvania  State  University. 
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4. 1.3.1  Initial  preparations. 

1.  Calibration  particles  -  In  order  to  control 
the  nuclei  concentrations  in  the  experiments, 
it  was  necessary  to  use  scientific  calibration 
particles.  To  model  the  cavitation  bubble 
nuclei,  latex  microspheres  [13]  were  used. 

These  came  in  the  form  of  a  suspension,  with 
approximately  10^  spheres  to  a  10-m£  bottle  if 
the  sphere  size  was  greater  than  24  microns  (pm) 
in  diameter.  The  number  was  10  if  the  sphere 
size  was  less  than  or  equal  to  24  pm.  To  model 
the  non-symmetr ic  solid  particulate  nuclei,  a 
calibrated  dust  particle  sample  was  used  which 
had  a  size  range  of  20  to  25  microns  in  mean 
diameter  guaranteed  by  analysis  [14], 

2.  PHA  calibration  curve  -  The  objective  here  was 
to  put  the  scattered  light  intensity  relation¬ 
ship  depicted  in  Figure  1  in  terms  of  actual  PM- 
tube  voltage.  This  was  done  by  using  an  alternate 
test  section  in  the  water  tunnel,  one  which  had 

an  open  top.  The  tunnel  was  subsequently  filled 
with  water,  and  latex  spheres  of  a  known  size 
were  then  manually  injected  into  the  laser  PV 
while  their  corresponding  voltages  were  recorded 
on  an  oscilloscope.  Only  a  few  sizes  needed 
to  be  tested  in  this  way  to  establish  the  entire 
voltage  calibration  curve. 
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3.  PHA  classification  channels  -  After  the  voltage 
calibration  curve  was  obtained,  it  was  then 
necessary  to  determine  the  size  and  number  of 
nuclei  classification  channels  in  the  PHA  that 
were  appropriate  to  the  experimental  conditions. 

An  attempt  was  made  to  measure  an  order  of  magni¬ 
tude  of  sizes  from  5  to  70  microns,  and  this 
range  was  spanned  with  10  classification  channels 
so  as  to  reduce  the  amount  of  experimentation 
necessary.  The  relative  widths  of  the  channels 
were  scaled  logarithmically;  that  is,  the  upper 
limit  voltage  amplitude,  A^,  of  a  given  channel 

i  was  related  to  the  upper  limit  voltage  amplitude 
of  the  next  channel,  A^+^,  according  to 

A.  .  -  A. 

i  =  8,  (4.1) 

i 

where  8  is  an  arbitrary  constant. 

4.  Matching  of  photo-multiplier  output  -  The 
matching  procedure  used  in  the  experiment 
followed  the  procedure  described  in  Section  4.1.2. 
Typically,  as  seen  in  Figure  10,  the  detector 
labeled  PM  1  was  the  one  used  to  establish  the 
voltage  calibration  curve,  and  then  detector 

PM  2  would  have  its  amplifier  gain  adjusted  until 
its  counting  time  was  matched  to  PM  1. 


Schematic  of  the  LSS  Equipment  Used  in  Experiments  to  Measure  Nuclei 
Distributions . 
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4. 1.3. 2  Typical  nuclei  test  procedure. 

1.  Tunnel  cleaning  and  preparation  -  The  tap 
water  source  was  connected  and  the  empty 
tunnel  was  allowed  to  fill  (refer  to  Figure  9). 
When  full,  the  drain  was  opened  and  the  pump  was 
turned  on  at  a  low  speed,  with  the  tap  water 
still  running.  The  tunnel  was  flushed  out  like 
this  for  about  20  minutes,  and  then  drained. 

Next,  the  pressure  line  was  connected  to 
the  distilled  water.  When  the  tunnel  was  full, 
the  pump  was  turned  on  and  the  water  was  allowed 
to  run  through  the  tunnel  for  about  20  minutes. 
After  this  water  was  drained,  the  tunnel  was 
refilled  with  a  fresh  sample  of  distilled  water. 

2.  Background  measurements  -  The  system’s  laser 
(spectra-Physics  4-watt,  argon-ion)  was  turned 
on  and  focused  into  the  test  section.  The 
laser's  blue  line  (4480  angstroms)  was  used, 
and  the  power  was  set  at  500  mW.  The  pump  was 
turned  on  and  brought  up  to  a  flow  of  12  feet 
per  second  at  the  test  section,  which  corresponds 
to  2000  rpm  at  the  motor.  The  background  noise 
level  remaining  in  the  tunnel  after  cleaning  was 
then  used  as  the  steady  count  rate  in  rechecking 
the  matching  of  the  photo -multiplier  signals. 
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When  rechecked,  both  PM  were  connected  to 
the  dual-input  combination  PHC/PHA  (refer  to 
Figure  10)  and  the  background  count  rate  was 
measured  10  consecutive  times.  This  was  done 
without  comparison  of  signals  in  the  PHC,  which 
is  an  effective  cutoff  rejection  ratio  of  R^  =  0.0 

3.  Nuclei  distribution  measurements  -  While  the 
tunnel  was  still  in  operation,  the  size 
concentration  to  be  measured  was  introduced  into 
the  tunnel  through  the  nuclei  injection  port 
(Figure  9)  .  Generally,  one  entire  bottle  of 
calibration  spheres  would  be  introduced  into 
the  tunnel  for  a  given  test.  After  the  spheres 
had  been  injected,  the  PHC/PHA  was  used  to  take 
10  consecutive  measurements  of  the  nuclei  count 
rate  with  no  comparison  in  the  PHC. 

4.  Begin  discrimination  effectiveness  data  - 

Without  stopping  the  tunnel,  the  PHC  was 

switched  to  the  discriminating  mode,  with  R£ 

set  at  0.50,  and  the  series  of  10  nuclei  counts 

was  taken.  This  procedure  was  repeated  for 

R  -  0.60. 
c 

5.  Holography  measurements  -  While  the  above 
measurements  were  being  taken,  the  KHS  was 
prepared  for  use.  After  the  last  nuclei  count 
was  made  with  R^  ■  0.60,  the  tunnel  pump  was 
turned  off,  and  two  or  three  holograms  were 
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immediately  taken.  Shutting  off  the  pump  was 
found  to  be  necessary  in  order  to  eliminate 
vibrations  that  blurred  the  hologram  image. 

Then  it  was  necessary  to  wait  approximately  one 
half  hour  for  the  holograms  to  be  developed 
before  proceeding  with  the  test. 

6.  Continue  discrimination  effectiveness  data  - 
After  the  holograms  were  developed,  the  tunnel 
was  restarted,  and  several  measurements  of  the 
count  rate,  using  no  comparison  in  the  PHC,  were 
made  as  a  recheck.  Then,  the  procedure  begun 

in  Step  4  was  resumed,  with  cutoff  rejection 
ratios  of  0.70,  0.80,  and  finally  0.90  being 
used  in  the  PHC. 

7.  Test  completion  -  When  the  measurement  of  the 
discrimination  effectiveness  data  was  completed, 
the  pressure  line  was  removed  from  the  distilled 
water  source  and  connected  directly  to  the  water 
tunnel  inlet,  and  the  tunnel  was  then  drained 
under  pressure. 

4. 1.3. 3  Scope  of  the  nuclei  tests.  The  typical  nuclei  test 
was  followed  for  a  number  of  different  nuclei  distributions.  To 
calibrate  the  area  matrix,  S^ ,  the  original  intention  was  to  measure 
calibration  sphere  monodispersions  representing  each  size  category  in 
the  PHA.  The  representative  monodispersons  were  chosen  to  be  the 
following:  5  pm,  7  pm,  10  ;a,  15.5  Pm,  19  pm,  24  pm,  31.2  va,  40  la. 
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50  Pm,  and  70  Dm,  representing  channels  1  to  10,  respectively. 

However,  during  the  process  of  conducting  the  experiments,  it  was 
found  that  the  holography  system  had  a  resolution  limit  of  15  pm; 
thus,  the  data  from  the  5-pm,  7-pm,  and  10-pm  tests  could  not  be 
interpreted. 

Two  mixtures  of  calibration  spheres  were  measured  using  this 
typical  test  procedure,  each  composed  of  three  different  sphere  sizes. 
Sphere  mixture  1  was  made  up  of  one  bottle  each  of  49.3-]jn,  31.2-pm,  and 
15. 5-pm  spheres;  sphere  mixture  II  was  composed  of  one  bottle  each  of 
50-pm,  24-pm,  and  15. 5-pm  spheres.  The  holograms,  taken  in  conjunction 
with  sphere  mixture  I,  could  not  be  interpreted,  so  only  the  data  from 
mixture  II  was  usable. 

Finally,  a  mixture  of  both  spheres  and  calibrated  dust  particles 
was  measured.  The  sphere  size  used  was  70  pm,  and  the  dust  was  the 
20  to  25  Pm  mean  diameter  particles  described  earlier. 

4. 1.3. 4  Hologram  reconstruction.  As  mentioned  in  Section 
4.1.2,  independent  measurements  of  the  concentrations  were  made  possible 
through  the  KHS.  As  can  be  seen  in  Figure  8,  such  measurements  were 
obtained  through  viewing  a  reconstructed  and  magnified  image  of  the 
hologram.  The  magnification  of  the  system  was  set  by  the  microscope 
objectives,  a  value  of  200x  for  this  investigation. 

As  the  hologram  was  moved  along  the  axis  of  the  reconstruction 
system,  different  planes  of  the  three-dimensional  image  would  come 
into  focus.  At  each  such  focus,  the  carriage  could  be  moved  horizontally 
or  vertically  to  allow  nuclei  to  be  counted.  Using  three  directions 
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of  motion,  the  number  of  spheres  in  an  arbitrary  volume  could  be 
counted  and  so  calculate  a  local  concentration. 

This  concentration  was  then  adjusted  by  a  correction  factor  that 
accounted  for  frequency-shifting  of  the  holography  laser  beam  in  water. 
This  correction  factor  was  calculated  by  Vikram  [15]  to  be 


actual 


(1.12)  V 


measured" 


(4.2) 


4.2  Experiments  Involving  Laser  Intensity  Distributions 

4.2.1  Objectives  of  Experiments.  After  the  experiments  involving 
laser  sphere  distributions  were  completed,  analysis  of  the  data  revealed 


did  not  satisfy  all  the  criteria  that  would  prove  the  validity  of  the 
inversion  scheme.  A  possible  cause  of  this  was  a  temporal  as  well  as 
spatial  variation  of  the  laser  intensity  distribution  at  the  PV. 

The  effect  of  this  would  be  that  each  individual  test  had  the  potential 
of  experiencing  a  different  matrix  of  area  elements,  even  though  the 
optical  arrangement  was  the  same  for  each  test. 

Experimental  measurements  of  the  laser  intensity  distribution  were 
obtained  so  as  to  accomplish  three  goals:  First,  to  document  spatial 
and  temporal  variations  in  the  laser  intensity;  second,  to  find  ways 
to  minimize  these  variations;  and  third,  to  isolate  the  contribution 
of  laser  intensity  non-uniformities  to  the  overall  inaccuracy  of  the 
count  rate  data. 

4.2.2  Description  of  Experiments 

4. 2. 2.1  Documenting  Intensity  variations.  As  shown  in 
Figure  11,  the  intensity  profile  was  measured  at  two  different  locations 
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along  the  laser  beam’s  path.  Position  1  was  2  to  3  feet  from  the  laser 
course  and  about  6  inches  in  front  of  the  pinhole  aperture.  Position  2 
was  approximately  6  feet  from  the  laser  source.  This  position  was 
supposed  to  be  a  model  of  the  actual  condition  of  the  PV  in  ARL's 
12-inch  diameter  water  tunnel;  thus,  Position  2  was  located  at  a  point 
where  the  beam  had  passed  through  a  1-inch  thick  plexiglas  window  and 
then  through  6  inches  of  water. 

Measurements  of  the  intensity  profile  were  obtained  through  the 
use  of  a  photo-detector  mounted  on  a  movable  carriage.  The  photo¬ 
detector  was  marked  with  a  30-vsn  diameter  pinhole  aperture,  and  the 
carriage  allowed  this  aperture  to  be  moved  horizontally  and  vertically 
in  fixed  increments  so  that  the  entire  cross-sectional  area  of  the  beam 
could  be  scanned . 

In  addition,  th^  ilrror  in  Figure  11a  could  be  removed  so  that 
the  beam  could  be  examined  on  a  screen  at  a  distance  of  approximately 
10  feet  from  the  laser  source.  The  intensity  profile  was  not  scanned 
with  the  photo-detector  at  this  distance.  Instead  a  lens  was  used  to 
magnify  the  beam  so  that  its  cross-section  could  be  inspected  visually 
on  the  screen. 

4. 2. 2. 2  Minimizing  variations.  When  the  mirror  was 
removed  and  the  beam  intensity  was  inspected  visually,  it  was  observed 
that  the  beam  became  seriously  distorted  within  18  inches  after  passing 
through  the  second  lens.  Therefore,  several  different  lenses  were 
used,  and  the  position  of  the  lens  with  respect  to  the  beam  was  varied 
In  order  to  improve  the  intensity  distribution. 


Figure  11a.  Schematic  of  the  Equipment  Used  in  Experiments  to  Measure  Laser  Intensity 
Distributions. 
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The  beam  was  observed  to  have  a  diffraction  pattern  around  its 
cross-section  which  was  rectangular  in  shape  like  the  laser  beam  itself. 
To  eliminate  this,  several  different  sizes  of  pinhole  apertures  were 
tried.  These  varied  from  25  Mm  in  diameter  to  400  Pm.  In  addition, 
a  series  of  rectangular  apertures  were  made  in  order  to  find  one  which 
would  most  closely  match  the  size  and  shape  of  the  beam's  central 
intensity  lobe  (see  Figure  lib). 

4. 2. 2. 3  Contribution  of  non-uniform  beam  to  count  rate 
inaccuracy.  It  was  decided  that  the  way  to  isolate  the  contribution  of 
the  PV  intensity  profile  to  the  inaccuracy  of  the  count  rate  distri¬ 
bution  given  by  the  PHA  was  to  obtain  measurements  of  a  monodispersion 
similar  to  the  tests  described  in  Section  4.1.3,  except  that  the  beam 
intensity  profile  would  be  measured  at  the  time  of  the  experiment. 

This  experiment  was  carried  out  using  70-pm  spheres  in  the  square 
tank  shown  In  Figure  11a.  The  water  in  the  tank  was  stirred  to  give 
the  spheres  a  velocity  through  the  control  volume.  Prior  to  the 
measuring  of  the  sphere  count  rate,  the  intensity  profile  was  measured 
in  both  Positions  1  and  2.  The  measurement  at  Position  2  was  made 
immediately  prior  to  measurement  of  the  count  rate. 

In  Appendix  C,  an  analysis  is  performed  which  generates  data 
allowing  a  comparison  to  be  made  between  the  count  rate  distribution 
obtained  from  this  experiment  and  the  distribution  that  should  have 
resulted  according  to  these  measurements  of  the  intensity  profile. 
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CHAPTER  5 

PRESENTATION  OF  DATA 


5 .1  Probe  Volume  Error  Analysis 


5.1.1  Experimental  Conditions.  In  Chapter  3,  the  formula  used 
to  estimate  errors  was  given  as 
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which  under  the  conditions  of  this  experiment  gave  rise  to  a  coincidence 
error  criterion  of 
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The  coefficient  in  Equation  (3.7)  is  the  combination  of  the  size  range 
parameter,  d  /d  .  ,  and  the  size  concentration  distribution  parameter, 
q.  The  size  range  programmed  into  the  PHA  was  5  pm  to  70  pm,  giving  a 
size  range  parameter  of  14.  The  value  of  q  was  obtained  from  a  composite 
plot  of  bubble  nuclei  concentration  measurements  compilea  by  Gates  [16]. 

A  linear  regression  of  the  data  on  his  log-log  plot  yields  a  functional 
relationship  approximated  as 

N(r)  -  Ar"2*25.  (5.1) 

Thus,  q  was  taken  to  be  equal  to  2.25,  and  when  substituted  into 
Equation  (3.3)  gave  the  relation  in  Equation  (3.7). 


The  lens  system  that  the  laser  beam  passed  through  produced  a 

width  of  0.67  mm  observed  at  the  Probe  Volume.  This,  together  with 

a  slot  aperture  width  of  0.025  inches  placed  over  the  PM-tubes,  gave 

-4  3 

a  PV  size  of  2.85  x  10  cm  . 


5.1.2  Analysis  of  Coincidence  Errors.  Given  the  size  of  the  PV 
used  in  the  experiments,  Equation  (3.7)  allows  the  probability  of 
coincidence  to  be  estimated  if  the  total  nuclei  concentration,  N^, 
is  known.  Table  1  shows  the  relationship  between  the  nuclei  concentration 
and  coincidence  probability  for  several  of  the  experiments  involving 
nuclei  distributions.  The  measurements  of  concentration  were  obtained 
by  using  the  holography  system. 

It  can  be  seen  that  the  tests  involving  the  40-pm,  50-pm,  and 
70-ym  spheres,  as  well  as  the  70-pm  spheres  and  dust  mixture,  all  had 
quite  low  coincidence  errors.  The  15.5-jjm,  19-pm,  and  sphere  mixture 
II  tests  all  had  very  high  coincidence  errors,  with  more  than  one  out 
of  every  three  events  being  a  coincidence. 

5.2  Performance  of  the  Dual-Detector  System 

5.2.1  Method  of  Proof.  As  detailed  in  Chapter  3,  the  proof  of 
the  dual -detector  system's  effectiveness  would  come  through  verifying 
the  theory  for  non-symmetric  particle  scattering  in  actual  practice. 

In  addition,  it  would  be  necessary  to  gain  an  estimate  of  the  asymmetry 
of  scattering  by  spheres. 

These  proofs  were  obtained  in  the  experiments  by  recording  the 
responses  of  several  calibration  sphere  monodispersions  to  different 
values  of  the  cutoff  asymmetry  ratio  (R^)  used  in  the  system's  Pulse- 
Height  Comparator  (PHC) .  In  the  same  way,  a  mixture  of  calibration 


TABLE  1 


RELATIONSHIP  BETWEEN  NUCLEI  CONCENTRATION  AND  COINCIDENCE  ERROR 


Nuclei  Test 

70-Vttn  spheres 
50-ym  spheres 
40-ym  spheres 
24-ym  spheres 
19-ym  spheres 
15.5-ym  spheres 

70-Un.  spheres  and  dust  particles 

sphere  mixture  II 

(50  Dm,  24  vtm,  and  15.5  ym) 


Total 

Concentration 

NT(in~3) 

Coinc idence 
Probability 

Pint 

407.7 

0.005 

835.7 

0.010 

3875 

0.047 

12510 

0.152 

30970 

0.377 

31310 

0.381 

3933.3 

0.048 

23896 

0.413 

V _ ! -  - -  -  — 


spheres  and  calibrated  dust  particles  was  tested  to  verify  the  theory 
for  non-symmetric  particle  scattering. 


5.2.2  Analysis  of  Data.  Table  3  gives  the  compiled  results  of 
all  the  experiments  involving  microsphere  distributions  and  their 
responses  to  changes  in  the  PHC  cutoff  rijection  ratio.  The  informa¬ 
tion  for  verifying  the  theory  for  non-symmetric  particle  scattering 
is  contained  in  the  data  from  the  70-pm  spheres  and  dust  mixture  test; 
however,  this  data  must  first  undergo  a  transformation  in  order  for 
it  to  be  in  terms  of  the  same  variables  used  in  the  theory. 

In  the  analysis  of  Kohler  and  Billet  [12],  the  amount  of  screening 
achieved  by  a  given  level  of  was  the  ratio  of  the  number  of  non- 
symmetric  (or  dust)  particles  rejected  to  the  total  number  of  dust 


In  this  investigation,  both  spheres  and  dust  particles  were  examined 
by  the  PHC,  so  that  the  rejection  ratio  given  was  actually  a  composite, 
defined  as 


N  .  .  .  +  N ,  .  , 

j  _  spheres  rejected  dust  rejected 

composite  N  ,  +  N, 

r  spheres  dust 


(5.3) 


In  the  tests  involving  only  microspheres,  the  rejection  ratio  of  the 
PHC  was  defined  in  yet  another  way  as 
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SPHERE  REJECTION 

Nuclei  Test 
70-ym  spheres 

50-um  spheres 

40-pm  spheres 

31.2-ym  spheres 

24-ym  spheres 

19-van  spheres 

15.5-ym  spheres 


TABLE  3 


DATA  AS  A  FUNCTION  OF  R  FOR  ALL  NUCLEI  TESTS 

c 


Cutoff  Asymmetry  Ratio 

In  PHC  (R  ) 
c 


PHC  Rejection  Ratio 


0.50 

0.747 

0.60 

0.818 

0.70 

0.870 

0.80 

0.911 

0.90 

0.952 

0.50 

0.533 

0.60 

0.638 

0.70 

0.804 

0.80 

0.912 

0.90 

0.919 

0.50 

0.751 

0.60 

0.895 

0.70 

0.926 

0.80 

0.957 

0.90 

0.968 

0.50 

0.690 

0.60 

0.  766 

0.70 

0.858 

0.80 

0.936 

0.90 

0.970 

0.50 

0.833 

0.60 

0.897 

0.70 

0.938 

0.80 

0.969 

0.90 

0.983 

0.50 

0.770 

0.60 

0.917 

0.70 

0.933 

0.80 

0.962 

0.90 

0.984 

0.50 

0.769 

0.60 

0.821 

0.70 

0.873 

0.80 

0.921 
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TABLE  3  (continued) 


SPHERE  REJECTION  DATA  AS  A  FUNCTION  OF  R  FOR  ALL  NUCLEI  TESTS 

c 


Cutoff  Asymmetry  Ratio 


Nuclei  Test 

in  PHC  (R  ) 
c 

PHC  Rejection  Ratio 

lO-ym  spheres 

0.50 

0.796 

0.60 

0.836 

0.70 

0.887 

0.80 

0.924 

0.90 

0.968 

70-pm  spheres  and 

0.50 

0.850 

dust  particles 

0.60 

0.879 

0.70 

0.912 

0.80 

0.947 

0.90 

0.975 

sphere  mixture  I 

0.50 

0.818 

(49.3  van,  36.2  van. 

0.60 

0.768 

and  15.5  pm) 

0.70 

0.925 

0.80 

0.919 

0.90 

0.960 

sphere  mixture  II 

0.50 

0.751 

(50  pm,  24  van,  and 

0.60 

0.885 

15.5  ym) 

0.70 

0.926 

0.80 

0.957 

0.90 

0.968 
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spheres 


^spheres  rejected 


N 


spheres 


(5.4) 


since  in  these  tests  there  were  no  dust  particles. 

Equation  (5.3)  can  be  written  in  terms  of  the  other  rejection 
ratio  as 
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or,  substituting  in  Equation  (5.2), 


composite 


R  .  +  (R,  .  M) 

sphere  dust _ 

1  +  M 


(5.6) 


where  M  is  defined  here  to  be  the  composition  ratio  of  the  mixture 
of  spheres  and  dust  particles. 


M  = 


N 

N 


dust 

spheres 


(5.7) 


Now,  expressing  the  rejection  ratio  in  terms  of  the  quantities  encountered 
in  this  investigation  gives 


dust 


<¥> 


composite 


A 

ID 


sphere" 


(5.8) 


The  composition  ratio  of  the  70-ym  sphere  and  dust  particle 


mixture  was  calculated  through  the  use  of  the  concentration  data  obtained 


* 
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from  holograms  during  that  test.  According  to  the  holography  measure¬ 
ments,  this  ratio  was 


2903.8  dust/in 

3 

1029.5  spheres/ in 


Using  this  value  of  M.  Table  2  gives  a  comparison  of  the  rejection 
ratio  in  Equation  (5.8),  as  functions  of  the  cutoff  asymmetry  ratio 
used  in  the  PRC. 


5.3  Performance  of  the  Matrix  Method 


5.3.1  Correlations  with  Sphere  Distributions 


5. 3. 1.1  Criteria  for  proof  of  matrix  method  validity. 

In  providing  the  validity  and  usefulness  of  the  matrix  method  as  a 
means  of  correcting  the  cavitation  nuclei  size  distribution  data 
given  by  the  light-scattering  system's  PHA,  two  proofs  are  required. 
First,  there  must  be  evidence  that  the  matrix  of  elements,  S^, 
representing  the  PV  are  intercepted  by  the  nuclei,  has  an  equi-diagonal 
character.  That  is,  the  experiments  must  give  equality  among  elements 
of  the  same  diagonal.  The  inversion  scheme  detailed  in  Reference  [10] 
predicts  that  such  a  behavior  will  occur  if  two  experimental  conditions 


are  met: 


1.  The  channel  widths  in  the  PHA  are  scaled 
logarithmically  (a  condition  detailed  in 
Section  4.1.3) . 

2.  The  area  function  is  independent  of  nuclei 
size  (see  Equation  (B.2)). 


i 


Second,  the  matrix  of  elements,  S  ,  obtained  in  experiments  should 
successfully  correct  the  PHA  data  from  an  experiment  involving  a  random 
size  distribution. 


5. 3. 1.2  Comparison  of  results  with  proof  criteria. 

Table  4a  shows  the  matrix  of  S  elements  that  were  obtained  through 
the  experiments  involving  nuclei  distributions.  Each  column  of  elements 
represents  the  results  of  one  test,  with  the  size  measured  listed  at 
the  head  of  each  column.  The  diagonals  are  shown  ordered  as 

—^io*  4b  lists  the  average  value  of  the  elements  in 

S £,  denoted  as  S^,  and  the  average  deviation  in  each  diagonal  from 
the  average  value.  From  those  two  tables  it  can  be  seen  that  the 
equality  of  diagonal  elements  does  appear,  and  that  the  range  of  spread 
in  the  deviation  from  each  diagonal  is  also  very  similar. 

With  the  area  matrix  exhibiting  the  characteristics  predicted 
by  the  inversion  scheme,  this  matrix  was  then  used  to  correct  PHA 
count  rate  data  from  two  different  random  nuclei  size  distributions. 

One  was  the  sphere  mixture  II,  and  the  other  was  an  actual  cavitation 
bubble  distribution  measured  in  the  12-inch  diameter  water  tunnel  at 
ARL.  In  each  case,  the  distribution  data  was  corrected  according  to 
Equation  (3.11)  for  three  different  matrices  of  elements.  One  was 
composed  of  the  actual  matrix  elements  obtained  by  experiments  and 
is  presented  in  Table  4a;  the  second  was  a  matrix  where  all  the  elements 
along  a  diagonal  were  equal  to  the  average  diagonal  value,  as  calculated 
in  Table  4b;  the  third  was  similar  in  form  to  the  second,  except  that 
the  value  along  a  diagonal  was  given  by  the  70-ym  sphere  test,  shown 
as  column  10  in  Table  4a. 


TABLE  4  a 


1.49x10 
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TABLE  4b 


SIMILARITY  OF  DIAGONAL  ELEMENTS  OF  S, . 

ij 


_  Average  Deviation  from 

Average  Value,  as  a  Z  of 


1.00  x  10 

1.37  x  10~4  4% 
9.08  x  10-5  122% 
3.02  x  10~5  99% 
1.83  x  10~5  58% 
1.31  x  10~5  78% 
6.28  x  10~6  62% 
3.68  x  10~6  57% 
2.00  x  10~6  67% 


1.17  x  10~6 


61% 


Tables  5a  and  5b  show  the  results  of  the  correlation  given  by 
each  of  these  matrices  compared  to  the  measurements  of  the  actual 
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concentration  obtained  by  holography.  None  of  the  matrices  show 
results  which  in  anyway  resemble  the  actual  concentrations.  The 
difference  in  matrix  elements  between  the  three  matrices  used  is  not 
large,  yet  the  difference  in  concentration  values  produced  is  very 

ft 

great,  suggesting  that  the  matrix-solving  routine  is  highly  sensitive 
to  errors  in  the  elements.  This  conclusion  is  strengthened  by  the 
observation  that  the  corrected  concentration  values  appear  to  get 
increasingly  inaccurate  as  one  goes  from  N^q  to  N^,  since  this  is  the 
order  in  which  the  concentrations  are  solved  for  in  the  inversion 
scheme. 


5.3.2  Correlation  with  Intensity  Distributions 

5. 3. 2.1  Comparison  of  intensity  distributions.  Since 
measurements  of  the  actual  PV  intensity  distribution  were  not  made 
during  the  experiments  involving  microsphere  distributions,  another 
series  of  experiments  was  conducted  to  investigate  this  phenomenon. 
Figures  12  through  14  show  the  results  of  laser  intensity  profile 
measurements  made  on  three  separate  occasions  without  changing  the 
laser  and  optical  arrangement.  In  each  figure,  part  A  shows  the 
intensity  profile  at  a  location  2  to  3  feet  from  the  laser  light  source, 
and  part  B  shows  it  at  about  6  feet  away  and  after  it  has  passed  through 
a  1-inch  thick  plexiglas  window  and  6  inches  of  water.  These  distri¬ 
butions  are  shown  normalized  with  respect  to  the  peak  intensity. 


TABLE  5a 


USING  Sy  MATRICES  TO  CORRECT  A  RANDOM  SIZE  DISTRIBUTION  OF 
NUCLEI  -  MIXTURE  SPHERE  SIZES 


Hologram 


Actual  Element  Column  Element  Average  Element 


0.0 

1.416 

X 

10' 

1.336 

X 

10' 

8.236 

X 

10 

0.0 

-1.641 

X 

io4 

-1.452 

X 

io4 

1.758 

X 

10 

0.0 

1.582 

X 

io4 

5.391 

X 

io3 

-1.318 

X 

10 

0.0 

-1.594 

X 

io3 

-2.869 

X 

io3 

-5.120 

X 

10 

1.600 

X 

io4 

-3.550 

X 

io3 

1.454 

X 

io3 

2.623 

X 

10 

0.0 

-3.980 

X 

io3 

-1.050 

X 

io3 

7.353 

X 

10 

1.468 

X 

io4 

4.903 

X 

io3 

9.973 

X 

io2 

-3.835 

X 

10 

0.0 

4.239 

X 

io2 

2.185 

X 

io2 

-3.640 

X 

10 

0.0 

5.050 

X 

io2 

2.414 

X 

io2 

1.061 

X 

10 

3.221 

X 

103 

2.315 

X 

io1 

2.905 

X 

io1 

-2.253 

X 

10 

0.0 

1.258 

X 

io1 

1.258 

X 

io1 

1.258 

X 

10 
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1H 


A 


TABLE  5b 


USING  S±j  MATRICES  TO  CORRECT  A  RANDOM  SIZE  DISTRIBUTION  OF 
NUCLEI  -  AIR  BUBBLES  IN  THE  ARL/PSU  12-INCH  WATER  TUNNEL 


Hologram 


Actual  Element  Column  Element  Average  Element 


0 

0.0 

-3.381 

X 

106 

2.273 

X 

io5 

7.713 

X 

10 

1 

1.441 

X 

IQ"1 

9.523 

X 

103 

1.191 

X 

io4 

1.371 

X 

10 

2 

7.205 

1.004 

X 

io4 

9.108 

X 

io3 

-1.056 

X 

10 

3 

3.602 

X 

101 

3.962 

X 

io3 

1.345 

X 

io3 

-5.854 

X 

10 

4 

6.484 

X 

io1 

1.349 

X 

io3 

3.152 

X 

io3 

2.041 

X 

10 

5 

2.161 

X 

101 

-3.696 

X 

io3 

-1.509 

X 

io3 

4.332 

X 

10 

6 

0.0 

1.114 

X 

104 

7.351 

X 

io2 

-3.015 

X 

10 

7 

0.0 

5.279 

X 

io2 

4.632 

X 

io2 

-2.855 

X 

10 

8 

0.0 

-4.604 

X 

io1 

7.185 

X 

io1 

8.312 

X 

10 

9 

0.0 

-5.010 

X 

io1 

-6.287 

X 

io1 

-1.769 

X 

10 

10 

7.205 

9.881 

X 

io1 

9.881 

X 

io1 

9.881 

X 

10 
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Figure  13a.  Normalized  Laser  Beam  Intensity  Distribution  Before 
Entering  Test  Section,  Obtained  on  10/7/81. 
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Figure  13b.  Normalized  Laser  Beam  Intensity  Distribution  Inside  Test  Section, 
Obtained  on  10/7/81. 
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Figure  14a.  Normalized  Laser  Beam  Intensity  Distribution 
Before  Entering  Test  Section,  Obtained  on 
10/16/81. 
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Figure  14b,  Normalized  Laser  Beam  Intensity  Distribution  Inside 
Test  Section,  Obtained  on  10/16/81. 
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It  can  be  seen  that  the  intensity  distributions  of  part  A  are  all 
very  similar  and  quite  uniform.  However,  the  same  cannot  be  said  for 
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their  counterparts  of  part  B.  Figure  15  brings  out  the  differences 
in  these  distributions  more  vividly  by  comparing  all  of  them  in  terms 
of  their  maximum  intensity  gradient  profiles.  Such  a  profile  illustrates 
the  largest  difference  in  illumination  that  a  scatterer  could  experience 
by  passing  through  any  given  portion  of  the  PV,  and  thus  is  an 
indication  of  the  PV's  non-uniformity.  The  intensity  gradient  is  the 
ratio  of  the  difference  in  intensity  across  an  interval  to  the  mean 
intensity  in  that  interval,  and  is  expressed  as  a  percentage  of  the 
mean  intensity  in  the  given  interval. 


distribution.  The  occurrence  of  these  differences  in  the  intensity 
distributions  suggested  the  possibility  that  optical  components  might 
be  having  an  effect.  It  was  observed  that  the  beam  began  to  depart 
from  a  uniform  intensity  distribution  as  the  beam  strayed  from  passing 
through  the  center  of  the  lens  system.  In  addition,  it  was  found 
that  the  diffraction  fringes  around  the  beam  could  be  eliminated 
through  the  use  of  a  square  aperture  that  was  50  pm  x  50  pm,  a  discovery 
made  through  successive  use  of  smaller  and  smaller  apertures  until 
the  diffraction  pattern  was  observed  to  disappear. 
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Comparison  of  Laser  Probe  Volume  Intensity  Distributions  in  Terms  of  Their  Maximum 
Intensity  Gradient  Profiles, 


CHAPTER  6 

DISCUSSION  OF  RESULTS 

6.1  Results  of  Coincidence  Error  Analysis 

A  very  likely  reason  for  the  differences  in  coincidence  error  in 

the  various  tests  is  the  difference  in  packaging  concentration  of 

calibration  spheres.  As  was  described  in  the  experimental  procedure, 

24  Um  marked  the  cutoff  size  where  calibration  spheres  larger  than 

this  came  10  ^  to  a  bottle,  and  sizes  smaller  than  this  came  10^  to  a 

bottle.  This  difference  in  the  order  of  magnitude  in  the  concentration 

is  reflected  in  the  difference  of  the  level  of  coincidence  error  because 

each  nuclei  test  used  an  entire  bottle  of  calibration  spheres. 

The  conclusion  from  these  observations  is  that  for  the  optical 

geometry  and  conditions  of  these  experiments,  the  highest  allowable 

nuclei  concentration  should  be  on  the  order  of  a  24-pm  sphere 

A 

concentration,  i.e.,  about  1.5  x  10  per  cubic  inch,  since  concentrations 
greater  than  this  result  in  coincidence  probabilities  greater  than  20%. 

In  the  analysis  of  the  nature  of  coincidence  errors  given  in  Chapter 
3,  it  was  mentioned  that  the  effect  of  nuclei  coincidences  in  the  PV 
could  cause  the  PM-tube  to  interpret  the  event  as  a  single  large 
nuclei,  so  that  one  would  expect  that  a  nuclei  distribution  with  a 
high  coincidence  probability  would  produce  a  PHA  count  rate  distribution 
with  a  disproportionate  number  of  counts  in  the  higher  channels. 

Table  1  Indicates  that  the  sphere  mixture  II  was  such  a  distribution, 
and  Table  5a  shows  that  the  concentration  calculated  by  matrix 
correction  of  the  count  rate  distribution  is  indeed  disproportionate. 
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predicting  scatterers  to  be  in  higher  size  categories  like  7,  8,  and 
10,  where  holography  shows  there  to  be  none.  Thus,  coincidence  errors 
could  contribute  to  the  inability  of  a  given  matrix  to  correct  a 
count  rate  histogram. 

6.2  Effectiveness  of  the  Dual-Detector  System 

6.2.1  Non-Symmetric  Particle  Rejection.  Kohler  and  Billet  [12} 
predict  that  if  a  concentration  of  non-symmetric  particles  whose  axial 
ratio  is  equal  to  1.67  are  analyzed  by  a  dual-detector  screening 
system  with  a  cutoff  asymmetry  ratio  set  at  R£  =  0.75,  84%  of  these 
particles  are  rejected  by  the  system.  If  one  interpolates  between 

R  =0.70  and  0.80  in  Table  2,  one  finds  that  the  actual  amount  of 
dust  particles  rejected  by  the  PHC  in  this  investigation  is  94%, 
indicating  that  the  theory  is  actually  conservative  in  its  estimate  of 
rejection  ratio. 

This  conclusion  is  strengthened  by  the  fact  that  the  dust 
particles  used  in  the  investigation  had  an  average  axial  ratio  of 
less  than  1.67,  around  1.3.  The  theory  predicts  that  as  the  axial 
ratio  decreases,  the  amount  of  rejection  should  also  decrease,  for  a 
given  value  of  R£  in  the  PHC.  Yet  in  actual  performance,  at  an  axial 
ratio  of  less  than  1.67,  the  amount  of  rejection  was  still  high. 

6.2.2  Sphere  Scattering  Asymmetry.  According  to  Table  3,  the 
amount  of  rejection  of  spheres  is  higher  than  expected.  As  the  cutoff 
rejection  ratio,  R£,  increases,  so  does  the  percentage  of  spheres 
rejected.  In  addition,  this  trend  appears  to  be  independent  of  size 
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and  concentration,  since  all  the  tests  show  nearly  the  same  amount  of 
rejection  for  a  given  value  of 

One  possible  reason  for  this  behavior  of  the  spheres  could  be 
the  influence  of  partial-illumination  errors.  In  Section  3.2.4,  it 
was  described  how  the  dual-detector  system  would  screen  out  spheres 
which  were  illuminated  non-unif ormly,  in  addition  to  screening  out 
non-uniformly  shaped  particles.  However,  the  relationship  between 
the  non-uniformity  of  a  sphere's  illumination  and  the  asymmetry  of 
its  scattered  light  intensity  was  not  known. 

Figure  15  gives  an  idea  of  the  range  of  illumination  non¬ 
uniformity  that  the  spheres  could  experience  in  the  PV.  As  an  example, 
consider  the  plot  corresponding  to  PV  measurement  of  10/16/81.  Here, 

72%  of  the  PV  width  has  a  maximum  illumination  non-uniformity  of 

greater  than  70%,  which  means  that  nuclei  passing  through  this  PV 

would  have  a  72%  chance  of  experiencing  greater  than  70%  non-uniform 

illumination  intensity.  If  it  is  assumed  that  70%  non-uniform  illumination 

of  a  sphere  results  in  70%  asymmetry  in  scattering  intensity,  then  a 

PHC  set  at  R^  =  0.70  would  reject  72%  of  the  nuclei  passing  through 

this  PV.  This  figure  is  on  the  same  order  as  the  percentage  of 

rejection  of  spheres  at  R  ■  0.70  in  Table  3. 

c 

Measurements  of  the  PV  intensity  distributions  experienced  by 
the  spheres  in  each  of  the  nuclei  tests  are  not  available.  Nevertheless, 
this  example  should  illustrate  the  fact  that  a  high  enough  illumination 
non-uniformity  can  exist  in  a  Probe  Volume  so  that  the  kind  of  sphere 
rejection  observed  in  Table  3  could  result. 


6.3  Effectiveness  of  the  Matrix  Method 


Through  observations  on  the  performance  of  the  matrix  method 
shown  in  Tables  5a  and  5b,  it  was  concluded  that  the  reason  for  the 
inability  of  the  given  matrices  to  correct  the  count  rate  data  must  be 
due  to  either  or  both  of  two  sources.  One  possibility  was  that  the 
actual  PV  intensity  distribution  (and,  therefore,  the  matrix) 
experienced  by  the  microspheres  was  different  in  each  test,  so  that 
the  Sy  matrix  measured  in  one  experiment  would  not  correct  count  rate 
data  from  a  different  experiment.  The  other  possibility  considered 
was  a  problem  of  counting  statistics;  that  is,  a  statistically 
insufficient  number  were  counted  in  the  tests.  Physically,  this  would 
mean  that  the  sample  was  not  large  enough  during  the  measurement  period  to 
give  a  true  picture  of  the  PV  intensity  distribution;  thus,  the 
resulting  S^.  matrix  could  not  give  a  complete  correction. 

The  information  in  Tables  6a  and  6b  gives  significant  insight  into 
the  influence  of  these  various  inaccuracies  on  the  performance  of  the 
matrix  method,  especially  the  comparison  of  corrected  count  rate  data 
from  sphere  mixture  II  given  in  Table  6b.  The  correction  based  on 
the  new  matrix  elements  represents  elements  which  have  been  modified 


to  account  for  the  influence  of  the  non-uniform  PV  illumination 


intensity.  That  the  correction  is  better,  yet  still  not  sufficiently 
accurate,  shows  that  the  non-unif ormity  of  the  PV  is  not  the  only  source 
of  inaccuracy. 

In  Appendix  C,  section  C.3,  a  description  is  given  of  a  phenomenon 
referred  to  as  the  "aperture  effect,"  which  is  shown  to  contribute 
to  the  inaccuracy  of  the  matrix  method  correction. 
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TABLE  6a 


COMPARISON  OF  COUNT  RATES  OBTAINED  FROM  EXPERIMENTS  WITH 
COUNT  RATES  OBTAINED  THROUGH  THE  INTENSITY  MAPPING  PROCEDURE 


Classification  Channel 
Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Count  Rate  from  Experi¬ 
ments  (nuclei/ second) 

342 

314 

159 

56 

35 

28 

27 

21 

15 

3 

Count  Rate  from  Intensity 
Mapping  (nuclei/sound) 

30 

20 

30 

30 

60 

60 

80 

100 

270 

320 

TABLE  6b 


COMPARISON  OF  S^  MATRIX  ELEMENTS  OBTAINED  FROM  EXPERIMENTS  WITH 
MATRIX  ELEMENTS  OBTAINED  THROUGH  THE  INTENSITY  MAPPING  PROCEDURE 


Column  10  Elements  Intensity  Mapping  Elements 


1.06 

X 

10-4 

1.19 

X 

io'5 

1.33 

X 

io-4 

7.95 

X 

IO’6 

4.63 

X 

IQ'5 

1.19 

X 

IO'5 

2.50 

X 

IO'5 

1.19 

X 

IO"5 

1.97 

X 

io"5 

2.38 

X 

IO"5 

1.42 

X 

io"5 

2.38 

X 

IO"5 

1.02 

X 

io'5 

3.18 

X 

IO"5 

5.84 

X 

io'6 

3.97 

X 

IO"5 

3.90 

X 

io-6 

1.07 

X 

IO'4 

1.49 

X 

io-6 

1.27 

X 

io"4 

CHAPTER  7 


CONCLUSIONS 


The  modifications  to  the  laser  light  scattering  system  developed 
by  Keller  have  been  shown  to  be  valid  and  relevant  means  of  solving 
inaccuracy  problems  in  this  system,  insofar  so  they  solve  the  inaccuracies 
they  were  designed  to  eliminate.  The  dual-detector  method  of  distin¬ 
guishing  microbubble  and  solid  particulate  cavitation  nuclei  has  been 
shown  to  exceed  expectations  for  screening  capability.  Yet  the  usability 
of  this  method  is  limited  by  a  different  problem,  that  of  spheres 
experiencing  significant  non-uniform  illumination  in  the  Probe  Volume. 

It  is  believed  that  this  is  the  primary  cause  of  such  a  large  per¬ 
centage  of  spheres  being  rejected  by  the  dual-detector  system  in 
addition  to  the  non-symmetr ic  particles. 

The  inversion  scheme,  or  matrix  method,  of  correcting  the  count 
rate  data  obtained  from  the  light-scattering  system's  Pulse-Height 
Analyzer  has  been  shown  to  be  relevant  to  the  problem  of  correcting 


for  the  distribution  of  light  intensity  in  the  PV,  since  the  matrix  of 
elements  obtained  from  experiments  conforms  to  the  predicted  equi- 
diagonal  form.  In  addition,  the  intensity-mapping  procedure  of 


Appendix  C  showed  that  true  values  of  the  correction  matrix,  S^,  could 
be  calculated  from  direct  measurements  of  the  PV  intensity  profile. 


This  investigation  also  revealed  another  possible  source  of 


inaccuracy,  that  of  insufficient  count  to  produce  a  statistically 
accurate  count  rate  distribution.  Although  the  extent  of  this  problem 
could  not  be  determined,  its  evidence  was  seen  in  that  the  inversion 
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scheme,  using  matrix  elements  generated  by  the  intensity-mapping 
procedure,  could  only  partially  correct  a  count  rate  distribution 
from  a  random  nuclei  size  concentration. 
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CHAPTER  8 

RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

Before  the  dual-detector  system  can  be  used  with  confidence  in 
actual  cavitation  tests,  the  problem  of  high  rejection  of  spheres 
must  be  solved.  An  investigation  should  be  conducted  into  the 
scattering  response  of  a  sphere  to  an  incident  beam  with  a  gradient 
of  light  intensity,  an  investigation  which  should  include  a  treatment 
of  the  Mie  scattering  equations  for  the  case  of  a  non-uniform  plane 
wave  incident  upon  a  spherical  particle.  The  results  of  such  an 
investigation  would  be  used  to  estimate  the  number  of  microbubble 
nuclei  that  would  be  rejected  by  a  dual-detector  system. 

In  terms  of  practical  additional  modifications,  a  mask  could  be 
placed  over  the  photo-multiplier  tube,  using  a  square  aperture 
instead  of  the  slit  aperture  currently  in  use.  Scans  of  the  PV 
intensity  distribution  would  reveal  the  region  with  the  least  non¬ 
uniformities,  and  the  size  and  position  of  the  aperture  could  be 
adjusted  to  examine  only  that  region. 

Concerning  the  problem  of  the  statistics  of  nuclei  counting,  an 
Investigation  should  be  conducted  as  to  how  many  nuclei  must  be 
measured  in  a  given  counting  period  to  insure  a  statistically  accurate 
sample.  Factors  like  the  nuclei  size  and  the  shape  of  the  FV  intensity 
profile  should  be  examined,  because  the  statistical  size  depends  on 
any  of  these  parameters. 
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CHAPTER  9 


SUMMARY 

This  investigation  has  presented  a  combination  of  modifications 
to  the  laser  light-scattering  system  designed  by  Keller  [2],  as  a 
means  of  solving  inaccuracy  problems  in  the  system  and  making  it  a 
reliable  means  of  measuring  nuclei.  The  two  primary  inaccuracy 
problems  were  the  inability  of  the  system  to  distinguish  between 
microbubble  and  solid  particulate  nuclei  and  the  persistent  registering 
of  nuclei  in  classification  channels  not  representing  their  true 
size. 

Theories  were  outlined  which  described  the  principles  governing 
these  inaccuracies,  and  modifications  to  the  Keller  system  were 
proposed  based  on  these  theories.  To  facilitate  the  screening  of 
solid  particulate,  or  non-bubble  nuclei  from  the  total  nuclei  distribution 
being  measured,  a  dual-detector  signal  analysis  system  was  proposed 
based  on  a  theory  of  scattering  by  non-symmetric  particles  in  the 
Mie  size  range.  To  correct  for  inaccurate  count  rates  reported  by 
the  system's  Pulse-Height  Analyzer,  a  mathematical  technique  called 
an  inversion  scheme  was  utilized.  A  formula  for  estimating  the  effect 
of  errors  due  to  nuclei  coincidence  was  presented,  based  on  the 
Poisson  probability  model. 

A  series  of  rigorous  experiments  was  designed  and  implemented  in 
order  to  test  the  effectiveness  of  these  modifications.  In  these 
experiments,  the  light -scattering  system  (LSS)  was  used  to  measure 
nuclei  distributions  in  a  highly  controlled  environment,  and  a 
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holography  system  was  used  to  calibrate  the  LSS  test  results.  In  this 
same  series  of  experiments,  microsphere  distributions  were  tested  to 
observe  their  response  to  various  levels  of  discrimination  in  the 
dual-detector  system. 

When  data  from  the  experiments  involving  nuclei  distribution 
failed  to  yield  a  matrix  of  elements  for  the  inversion  scheme  that 
could  successfully  correct  random  nuclei  distributions,  a  second 
series  of  experiments  was  conducted  to  investigate  changes  in  the 
intensity  profile  of  the  Probe  Volume. 

In  the  process  of  analyzing  the  combined  results  of  these 
experiments,  another  source  of  inaccuracy  in  the  light-scattering 
system  was  encountered,  which  was  the  problem  of  the  statistics  of 
nuclei  counting,  or  not  measuring  enough  nuclei  in  a  given  counting 
period  to  give  a  statistically  accurate  count  rate  distribution. 

Evidence  of  the  problem  was  discovered  through  the  use  of  the  intensity¬ 
mapping  technique,  which  was  introduced  in  this  investigation  as  a 


means  of  determining  the  elements  of  the  area  matrix,  S  ,  by  direct 
measurement  of  the  PV  intensity  profile.  However,  the  means  for 
eliminating  this  problem  were  not  developed  in  this  investigation. 

The  conclusion  of  the  investigation  was  that  both  modifications 


to  the  Keller  system  were  relevant  to  the  inaccuracy  problems 


observed,  and  that  their  capability  of  solving  these  inaccuracies  has 
been  proven,  but  that  their  usability  in  actual  cavitation  experiments 
would  be  limited  until  the  problem  of  counting  statistics  and  non¬ 


uniform  illumination  of  spheres  in  the  PV  were  solved.  Recommendations 


were  made  for  investigating  these  problems. 


-l 


REFERENCES 


1.  Holl,  J.  W.,  "Nuclei  and  Cavitation,"  ASME  Transactions,  J. 

Basic  Engineering,  December  1970,  pp.  681-688. 

2.  Morgan,  W.  B.,  "Air  Content  and  Nuclei  Measurements,"  Report  of 
the  13th  ITTC  Cavitation  Committee,  1972. 

3.  Keller,  A.  P.,  "The  Influence  of  the  Cavitation  Nucleus  on 
Cavitation  Inception,  Investigated  with  a  Scattered  Light  Counting 
Method,"  ASME  Transactions,  J.  Basic  Engineering,  December  1972. 

4.  Billet,  M.  L.,  and  Gates,  E.  M. ,  "A  Comparison  of  Two  Optical 
Techniques  for  Measuring  Cavitation  Nuclei,"  International 
Symposium  on  Cavitation  Inception,  ASME,  New  York,  New  York, 
December  1979. 

5.  Treaster,  A.  L. ,  "Cavitation  Hysteresis,"  M.S.  Thesis,  The 
Pennsylvania  State  University,  June  1964,  pp.  12-15. 

6.  Gates,  E.  M.,  et  al.,  "Cavitation  Inception  and  Nuclei  Distri¬ 
butions  -  Joint  ARL/CIT  Experiments,"  CIT  Technical  Report 

No.  E244.1,  September  1979,  p.  23. 

7.  Yungkurth,  C.  B.,  "Optimization  of  a  Light -scatter ing  System,” 

B.S.  Thesis,  The  Pennsylvania  State  University,  May  1979. 

8.  Avellan,  F.,  and  Resch,  F. ,  "A  Scattering  Light  Probe  for  the 
Measurement  of  Oceanic  Air  Bubble  Sizes,"  International  Journal 
of  Multi-phase  Flow,  1981. 

9.  Ling,  S.  C.  and  Gowing,  S. ,  "Measurements  of  Microbubbles  in 

a  Water  Tunnel,"  19th  American  Towing  Tank  Conference,  University 
of  Michigan,  Ann  Arbor,  Michigan,  July  1980. 

10.  Holve,  D.,  and  Self,  S.  A.,  "An  Optical  Particle-sizing  Counter 
for  In-situ  Measurements,"  Technical  Report  for  Project  SQUID, 
Subcontract  8960-7,  January  1978. 

11.  Keller,  A.  P. ,  et  al.,  "Comparative  Investigations  of  the 
Scattered  Light  Counting  Method  for  the  Registration  of  Cavitation 
Nuclei  and  the  Coulter  Counter,"  Prepared  for  National  Science 
Foundation,  Michigan  University,  January  1974. 

12.  Kohler,  R.  A.,  and  Billet,  M.  L. ,  "Light  Scattering  by  a  Non- 
spherical  Particle,"  Cavitation  and  Polyphase  Flow  Forum  -  1981, 
ASME,  Boulder,  Colorado,  June  1981. 

13.  (Reference  to  Duke  Scientific  microspheres  -  p.  20)  Duke 
Scientific  Corporation,  445  Sherman  Avenue,  Palo  Alto,  California, 
94306. 


I 


81 


14.  Civil  Engineering  Department,  Oklahoma  State  University, 

Stillwater,  Oklahoma,  74047. 

15.  Vikram,  Chandra  S.  ,  Research  Associate,  Materials  Research 
Laboratory,  The  Pennsylvania  State  University,  private  communi¬ 
cations  . 

16.  Gates,  E.  M. ,  and  Bacon,  J. ,  "A  Note  on  the  Determination  of 
Cavitation  Nuclei  Distributions  by  Holography,"  SNAME  J.  Ship 
Research,  Vol.  22,  No.  1,  March  1978 

17.  Van  de  Hulst,  H.  C. ,  Light  Scattering  by  Small  Particles,  New  York, 

Wiley,  1957.  ~ .  *  "  ' 

18.  Kerker,  M. ,  Scattering  of  Light .  New  York,  Academic  Press,  1969. 

19.  Billet,  M.  L.  Private  Communications. 


I 


''AD- A 1 27  946  A  L IGHT -SCATTERING  SYSTEM  TO  MEASURE  CAVITATION  NUCLEI 
ANALYSIS  AND  CAll-.(U)  PENNSYLVANIA  STATE  UNIV 
UNIVERSITY  PARK  APPLIED  RESEARCH  LAB..  C  B  YUNGKURTH 
UNCLASSIFIED  30  SEP  82  ARL/PSU/TM-82-204  F/G  20/6 


APPENDIX  A 


THEORETICAL  DEVELOPMENT  OF  THE  OPTIMUM  MEASURING  OPTICS  LOCATION 


The  basis  for  a  theoretical  description  of  this  light-scattering 
phenomenon  comes  primarily  from  the  work  of  G.  Mie  in  1909  on  the 
scattering  of  coherent  light  from  spherical  particles,  as  derived  from 
Maxwell’s  equations  of  the  electromagnetic  character  of  light  [17,18]. 
A  brief  summary  of  this  development  follows. 


A.l  Mie  Theory  Derivation 

The  Mie  theory  of  light  scattering  from  a  spherical  particle 
arises  from  the  solution  of  the  scalar  form  of  Maxwell's  wave  equation 
in  spherical  polar  coordinates: 


1 

32(P") 

+  1  3 

-ln0  Hi 

p 

_3p2  _ 

p^iae  w  1 

sin0  3Q 

p2sin20  9(f> 


a2-*  2 

— 2  +  k  11 


0, 


(A.l) 


where  7r  is  a  scalar  component  of  a  wave  vector  defined  by  the  electric 
*►  -► 

vector  E,  and  the  magnetic  vector,  H.  Equation  (2.1)  is  separable  into 
three  functions  as 


R(p)  0  (0)  4>  (4>) 


(A.  2) 


The  resulting  equations  are  familiar  ordinary  differential 
equations  which  have  the  following  solutions: 
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t  : 


0(6)  =  P  vra;  (cos6) 
n 


♦  (0)  =  sin(m<{>),  cos(m<(>) 


(A.  A) 


(A.  5) 


where  J  and  N  are  the  half-integral  order  Bessel  functions  and  P 

n 

are  the  Legendre  polynomials,  with  the  restriction  that  -n  _<  m  <  +n. 

The  general  solution  for  the  scalar  wave  equation  is  then  the  linear 
superposition  of  all  the  particular  solutions: 


plT  =  Jo  m=-n  [CnVkP)  +  dn*n(kP)  1  (cOS0)  J  [amCOs(me> 


+  b  sin(m<J>)  ] 

m 


(A.  6) 


In  regards  to  the  problem  of  scattering  from  a  sphere,  two  cases 
must  be  represented:  The  equation  of  the  incident  wave  and  that  of 
the  scattered  wave.  A  summation  of  the  above  expression  over  m  gives 
these  equations: 


p7r  1  =  -L  E  in_1  ((2n+l)/n(n+l)J  (kp)  P  (1)  (cos0)cos<f>  (A.7) 
1  kZ  n-1  n  n 


incident  wave 
components 
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p*,2  *-4  j  ln-1  l(2n+l)/n(n+l)]  a  C  (kp)P  (k)  (cos0)cos<f>  (A.9) 
1  kZ  n-1  n  n  n 


scattered  wave 
components 


CO 

pir*  *  -  Z  i”"1  t(2n+l)/n(n+l)]  b  T,  (kp)P  (1)  (cos9)  sin(J>.  (A. 10) 
2  k2  n-1  n  n  n 


C^(kp)  =  ^(kp)  +  i^Ckp)  appears  in  the  scattered  wave  expression 

because  it  vanishes  at  infinity,  a  necessary  property  of  the  scattered 

wave.  In  addition,  the  coefficients  a  and  b  (not  to  be  confused 
*  n  n 

with  a  and  b  )  appear  as  the  arbitrary  factors  which  characterize 
n  cn 

the  scattered  wave.  They  are: 


V<y>*n(x)  "  ®^nCy)^n'(x) 
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m^n'(y)^n(x)  -  ^n(y)^n'(x) 
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(A. 12) 


relative  refractive  index 
where  x  *  size  parameter  =  — ~  y  =  mx. 

Because  the  function  £(kp)  vanishes  at  infinity,  it  can  be  approximated 

-ikp* 

by  a  function  of  the  type  e  .  Thus,  one  final  simplifying  change 

is  made  in  the  expressions  for  the  scattered  wave: 


pir.8  *  (^)  e  ikp  I  [(2n+l)/n(n+l)2)  [a  Tr  (cos0) 
1  kp  n-1  n  n 


+  b  T  (cos8)  ]  coa<P 
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^Assuming  that  p  is  always  measured  in  the  far  field. 
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and 


P*  2  =  (■=£)  e_ikp  I  [(2n+l)/n(n+l)  ]  [b  tt  (cosS) 

s  kp  _ .  n  n 

n=i 
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where,  for  the  purposes  of  calculation,  tt  and  T  are  defined  as 

n  n 


p  (cos9)  =  P  (cos6) /sin9 
n  n 


(A. 15) 


T  (cos9)  -=  A  [P  (1)(cos9)l 
n  don 


(A. 16) 


Now,  the  scattered  electric  field  vector,  and  subsequently  the 
scattered  light  intensity,  tnay  be  defined. 

The  general  form  of  the  vector  wave  equation  is 


iwt 


(A. 17) 


Having  obtained  the  scalar  factor  A  from  the  solution  of  the  scalar 
wave  equation,  the  scattered  components  of  the  electric  field  can  be 
expressed  as: 


and 


E9  *  (^>  e1(wt_kp)  cos*  S2(9) 
E9  -  (^)  ei(wt_kp)  sin*  S1(9) , 


(A. 18) 


(A. 19) 


where  S^(9)  and  S2(9)  are  known  as  the  complex  amplitudes  of  the 


scattered  wave,  and  are  defined  as 
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Therefore,  substituting  Equations  (A. 25)  and  (A. 26)  into  Equation  (A. 22) 
gives  the  final  expression  for  the  intensity  as 


Itotai  =  ~ ~2~~2  l<sin<kp)  _  icos(kp))  sin^S^O)  -  (sin(kp+6) 
k  p 

-  icos(k+6))  cos^CS)  .  (A.  28) 

A. 2  Characteristics  of  Mie-Scattered  Light 

It  can  be  seen  from  section  A.l  that,  in  general,  this  scattered 
light  intensity  can  be  summarized  in  function  form  as 

I/IQ  -  f(d,  m,  p,  9,  <{>,  X),  (A. 29) 

where  I  **  scattered  light  intensity 

IQ  =  incident  light  intensity 
d  »  spherical  particle  diameter 

m  *  relative  refractive  index  of  particle  with  respect 
to  medium 

p  *  distance  from  particle  location  to  intensity 
measurement  location 

6,$  *  solid  angle  into  which  light  is  scattered 
X  ■  laser  light  wavelength. 

This  scattered  light  intensity  shows  a  strong  angular  dependence,  and 
it  is  because  of  this  dependence  that  an  optimum  position  (6,<|>)  of 
the  receiving  optics  was  investigated. 


APPENDIX  B 


SUMMARY  OF  THE  INVERSION  TECHNIQUE 

B.l  Introduction 

In  the  process  of  developing  the  ARL/PSU  Laser  Light-scattering 
System  (LSS)  for  the  measurement  of  microbubble  and  solid  particulate 
in  water.  Reference  17]  shows  that  such  nuclei  will  scatter  light 
in  the  manner  predicted  by  the  calculations  of  Mie  in  1909.  This  was 
done  by  observing  the  response  function  of  single  particles  as  they 
passed  through  a  beam  of  laser  light.  However,  experience  has 
consistently  shown  that  when  this  same  experimental  technique  is  used 
to  measure  concentrations  of  nuclei  in  a  continuous-flow  situation, 
the  technique  fails  to  give  an  accurate  representation  of  either  the 
nuclei  concentration  or  the  nuclei  size  distribution. 

In  January  of  1978,  Holve  and  Self  [10]  were  developing  a  similar 
LSS  and  encountered  the  same  experimental  difficulties.  Their 
solution  was  to  develop  a  numerical  method  of  evaluating  the  nuclei 
counter  data,  based  on  a  correction  for  the  particle  trajectory 
through  different  parts  of  the  LSS  measuring  control  volume.  Mathe¬ 
matically,  this  correction  can  be  described  as  follows. 

Previous  to  this  numerical  scheme,  the  amplitude  of  the  voltage 
pulse  corresponding  to  a  particle  entering  the  measuring  control 
volume  (CV)  was  thought  to  be  given  by 


A  -  C  IQ  F(o  ,  fl) 


(B 
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where  A  =  voltage  amplitude 

G  *  gain  of  the  photo -multiplier  and  associated 
electronics 

Iq  =  intensity  of  the  laser  light  in  the  CV 
F(a,£2)  =  scattering  response  function 
a  -  particle  size  parameter 

3  solid  angle  into  which  light  is  scattered. 

According  to  Holve  and  Self,  Equation  (B.l)  should  be  modified  to  read 

A  =  G  I0(y,z>  H(yfz,ct)  F(ot,ft),  (B.2) 

where  H  «  dimensionless  function,  varying  between  0  and  1, 

describing  the  percentage  of  scattered  light 
from  the  CV  that  reaches  the  PM, 

and  where  I  is  now  seen  to  be  a  function  of  the  coordinates  transverse 
o 

to  the  flow  direction  (see  Figure  3).  It  can  be  seen  that  Equation  (B.2) 
will  collapse  to  Equation  (B.l)  if  the  CV  illumination  intensity  (Iq) 
is  uniform,  and  if  the  transfer  function  (H)  is  equal  to  unity. 

Similarly,  it  was  previously  known  that  the  nuclei  counting  rate 
is  related  to  the  nuclei  number  density  by 

Ci  "  USNi»  (B,3) 

where  C^  *  nuclei  counting  rate 

U  ■  mean  flow  velocity 
S  »  area  of  CV  normal  to  U 


■  nuclei  number  density. 


rr  i 


Holve  and  Self  modified  Equation  (B.3)  to  provide  the  matrix  equation 


C,  -  US. .  N. , 
i  i]  J 


(B.4) 


where  *  C(A^)  =  column  matrix  of  signal  count  rates  for 

signal  amplitudes  in  the  range  to  a^  +  AA^ 

Nj  =  N(F^)  =  column  matrix  of  nuclei  number  densities  for 
nuclei  in  the  size  parameter  range  that  yield 
normalized  response  functions  in  the  range  F^  to 
F  +  AF 

j  3 

S  *  SCA^,  Fj)  *  equivalent  cross-sectional  area  of  the 
PV  which  yields  normalized  amplitudes  in  the  range 
A^  to  A^  +  AA^  for  particles  having  normalized 
response  functions  in  the  range  to  F^  +  AF^ 
i,j  =  indices  referring  to  the  channel  number  on  the  pulse- 
height  analyzer  (PHA) . 

The  actual  data  received  in  a  nuclei-counting  experiment  are  the 
collected  in  the  channels  of  the  PHA.  The  information  desired  is 
the  nuclei  number  density,  ,  in  each  size  range,  j,  of  the  PHA.  Thus, 
Equation  (B.4)  is  rewritten  as 


N  «  —  C  S-1 
j  U  i  ij  ’ 


(B.5) 


1  , 
'  /*"! 


where  *  inverse  of  the  matrix  in  Equation  (1.4). 

Thus,  it  can  be  seen  from  Equation  (B.5)  that  the  nuclei  number  density 


in  one  channel  is  actually  a  weighted  sum  of  the  count  rates  from 
every  channel.  This  is  the  basic  proposition  of  the  numerical  method 


of  Holve  and  Self. 
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B.2  Determination  of  the  Area  Matrix 

It  is  obvious  from  Section  B.l  that  the  key  to  determining  the 
nuclei  number  densities  is  the  ability  to  determine  the  area  matrix, 
S^.  Recalling  Equation  (B.2),  it  can  be  seen  that  this  area  is 
related  to  the  beam  intensity  and  transfer  function,  and  therefore  it 
can  have  a  dependence  on  particle  size  in  the  most  general  case. 

Therefore,  the  most  direct  and  the  most  certain  method  of 
determining  the  area  matrix  would  proceed  as  follows.  A  monodisperse 
particle  distribution  of  known  size,  d^,  and  concentration,  N(dk> , 
is  passed  through  the  measurement  volume  at  a  known  speed,  U.  Signal 
pulses  are  accumulated  in  the  PHA  for  a  specified  length  of  time, 
yielding  a  count  rate  distribution,  The  matrix  Equation  (B.4) 

then  reduces  to 


C  -  U  S,.  N.  , 
i  ik  k’ 


(B.6) 


since  N(dj)  ■  0  for  all  j,  except  j  =  k.  Thus,  Equation  (B.6)  will 
give  all  the  i  elements  of  column  j  **  k  of  the  matrix.  The  procedure 
is  then  repeated  for  the  j  size-range  categories  of  the  PHA. 

This,  however,  can  be  a  time-consuming  project.  Holve  and  Self 
have  shown  that  this  procedure  can  be  greatly  simplified  if  two 
criteria  are  met: 

(1)  The  amplitude  categories  in  the  PHA  are  scaled  logarithmi¬ 
cally. 

(2)  The  area  function,  specifically  H(x,z,a),  is  independent 
of  size  parameter  a. 


*1 
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t  : 


However,  given  that  the  necessary  criteria  are  met,  the  set  of 
equations  in  Equation  (B.6)  reduced  further  to 


Ci  "  U  S*  \  <B‘7> 


S2,  “  UN  Ci’  (B.8) 

m 


J 


N  *  concentration  of  a  particle  size,  d  ,  from  the 
m  r  max ’ 


where 

m 

largest  size  class  in  the  PHA. 

The  implication  of  Equation  (B.8)  is  twofold.  First,  it  implies 
that  the  form  of  the  S  matrix  is  such  that  all  the  elements  on  any 
diagonal  are  equal.  Second,  it  implies  that  all  the  elements  of  the 
matrix,  i.e.,  H  =  1,  m,  can  be  determined  from  a  single  calibration 

experiment,  using  a  concentration  of  a  particle  size  d  . 

m 

Once  the  elements  of  the  S  matrix  have  been  determined,  the  only 
step  remaining  is  to  calculate  its  inverse,  Sy .  This  should  be  a 
relatively  simple  matter,  since  will  always  be  in  upper -triangular 
form  with  non-zero  elements. 

For  a  given  geometry  and  flow  conditions  of  an  experimental 
setup,  this  matrix  need  only  be  determined  once.  Thereafter,  the  nuclei 
concentrations  may  be  calculated  directly  from  the  PHA  count-rate  data, 
according  to  Equation  (B.5) . 
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APPENDIX  C 

THE  INTENSITY  MAPPING  PROCEDURE 

C.l  Description  of  the  Analysis 

In  Section  4. 2.2. 3,  it  was  mentioned  that  measurements  of  a 
monodispersion  of  70-pm  spheres  were  made  in  conjunction  with  measure¬ 
ments  of  the  Probe  Volume  (PV)  intensity  profile.  These  data  were 
then  used  to  make  a  comparison  between  the  count  rate  distribution 
that  was  measured  experimentally  and  the  distribution  that  should  have 
resulted  according  to  the  intensity  profile  that  was  measured 
experimentally. 

Figure  16  demonstrates  how  an  analysis  referred  to  as  the 
intensity  mapping  procedure  was  used  to  accomplish  this  comparison. 

This  analysis  assumes  that  a  spatially  uniform  nuclei  concentration 
of  only  one  nuclei  size  (i.e. ,  a  monodispersion)  is  being  measured 
and  that  all  nuclei  have  an  equal  probability  of  sampling  all  portions 
of  the  PV. 

Figure  16  is  expressed  in  terms  of  normalized  quantities.  The  right- 
side  plot  is  the  light-scattering  calibration  curve  of  Figure  1. 

Instead  of  nuclei  size,  the  right-side  abscissa  is  a  fractional  nuclei 
size,  where  the  normalized  nuclei  size  is  expressed  as  a  percentage 
of  the  monodispersion  nuclei  size.  Likewise,  the  right-side  ordinate 
is  a  normalized  scattered  light  intensity.  Considering  the  left-side 
plot,  this  is  an  intensity  profile  obtained  by  measurements  of  the 
laser  light  intensity  at  the  center  of  the  PV.  The  measurements  were 
in  a  direction  transverse  to  the  direction  of  the  incoming  nuclei,  a 
direction  which  is  parallel  to  the  beam's  width.  Thus  the  left-side 
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abscissa  is  a  normalized  PV  width.  Since  the  assumption  was  made  that 
all  nuclei  sample  the  PV  equally,  the  abscissa  takes  on  an  additional 
meaning;  namely,  that  a  fraction  of  this  PV  width  is  equal  to  a  fractional 
probability  of  nuclei  sampling  that  portion  of  the  PV  width.  Again, 
the  ordinate  of  this  curve  is  likewise  a  normalized  intensity,  where 
here  the  normalizing  quantity  is  the  peak  laser  light  intensity  in 
the  PV. 

Having  defined  the  quantities  in  Figure  16,  the  intensity-mapping 
analysis  proceeds  as  follows: 

1.  The  nuclei  classification  channel  divisions  (as  programmed 
into  the  PHA)  are  normalized  by  the  monodispersion  nuclei 
size.  In  this  particular  case,  it  is  assumed  that  the 
experiment  had  a  monodispersion  of  70-pm  spheres.  The  70-um 
size  happens  to  be  from  the  highest  size  category  in  the  PHA, 
so  that  when  the  channel  divisions  are  normalized,  all  10 
channel  divisions  will  appear.  Had  the  monodispersion 
corresponded  to,  for  instance,  the  eighth  size  category,  then 
when  the  channel  divisions  were  normalized  only  channels 

1  through  8  would  have  appeared  on  the  right-side  abscissa 
of  Figure  16. 

2.  The  calibration  curve  (right-side  plot)  relates  a  set  of 
normalized  nuclei  size  ranges  to  a  set  of  normalized 
scattered  light  intensity  ranges. 

3.  The  set  of  normalized  scattered  light  intensity  ranges 
determine  a  set  of  normalized  PV  intensity  ranges.  There 
is  an  assumption  here  which  can  be  expressed  thusly:  when 
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a  nucleus  is  illuminated  by  some  fraction  of  the  PV  peak 
intensity,  the  optical  recording  system  "sees"  this  as  a 
nucleus  of  some  fraction  of  its  true  size  being  illuminated 
by  the  PV  peak  intensity. 

4.  The  PV  intensity  profile  relates  a  set  of  normalized  PV 
intensity  ranges  to  a  set  of  fractions  of  the  PV  width. 

5.  The  set  of  fractions  of  the  PV  width  is  equal  to  the  set 

of  probabilities  of  nuclei  experiencing  the  set  of  normalized 
PV  intensity  ranges. 

6.  The  set  of  probabilities  gives  the  set  of  numbers  of  nuclei 
which  experience  the  set  of  normalized  PV  intensity  ranges. 

C.2  Application  of  the  Analysis 

When  the  various  physical  quantities  involved  are  considered  in 
the  way  outlined  by  this  intensity  mapping  analysis,  the  intensity¬ 
mapping  procedure  can  be  carried  out  in  the  order  given  above,  or  in 
a  "reverse"  order.  In  the  order  stated  above,  the  analysis  gives  the 
number  of  nuclei  which  will  appear  in  each  classification  channel  out 
of  a  given  sample  size  (like  1000  nuclei),  for  a  known  PV  intensity 
distribution. 

Table  6a  gives  the  results  of  this  analysis,  compared  to  the 
actual  distribution  of  counts  in  the  classification  channels  that 
occurred  for  this  experiment.  In  Table  6b  the  count  rates  obtained 
by  this  analysis  are  taken  a  step  further  through  calculation  of  a 
set  of  equivalent  matrix  elements,  S^,  based  on  the  experimental 
conditions  that  yielded  the  column  10  elements  of  Table  5.  Finally, 
in  Table  6c  these  data  are  taken  yet  another  step  further  by  using  the 
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TABLE  6c 

COMPARISON  OF  THE  CORRECTION  GIVEN  BY  EXPERIMENT  MATRIX  ELEMENTS 
WITH  THE  CORRECTION  GIVEN  BY  INTENSITY  MAPPING  ELEMENTS 

3 

Corrected  Sphere  Mixture  II  (nuclei/ in  ) 


Concentration 

From  Intensity 
Mapping  Elements 

From  Column  10 
Elements 

Holography  Concen¬ 
tration  Measurements 

N1 

0.15 

1.452  x  10A 

0.0 

N2 

0.60 

5.391  x  103 

0.0 

N3 

3.75 

-2.869  x  103 

0.0 

N4 

8.94 

1.454  x  103 

1.600  x  10A 

N5 

24.37 

-1.056  x  103 

0.0 

N6 

28.49 

9.973  x  102 

1.468  x  10A 

N7 

47.36 

2.185  x  102 

0.0 

N8 

56.81 

2.414  x  102 

0.0 

N9 

96.13 

2.905  x  101 

3.221  x  103 

N10 

55.58 

1.258  x  101 

0.0 
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matrix  elements  of  Table  6b  to  calculate  a  corrected  sphere  concentration 
distribution.  Based  on  the  count  rate  data  obtained  from  the  sphere 
mixture  II  experiment.  It  can  be  seen  that  the  concentrations 
calculated  with  these  equivalent  matrix  elements  are  more  reasonable 
than  those  obtained  with  the  matrix  elements  derived  from  the  70-pm 
sphere  count  rate  measurements  (i.e.,  the  column  10  elements). 
Nevertheless,  they  still  are  not  sufficiently  accurate,  as  is  seen 
when  they  are  compared  to  the  actual  concentration  data  obtained 
through  holographic  measurements. 

The  procedure  could  also  be  performed  in  a  "reverse"  order;  that 
is,  an  analysis  which  would  show  the  PV  intensity  distribution  that 
would  produce  a  known  distribution  of  nuclei  in  the  PHA  classification 
channels.  Such  an  analysis  was  performed,  and  the  result  is  the  curve 
that  appears  to  be  approaching  a  delta-function  on  the  left-side  plot 
of  Figure  16. 

C.3  The  Aperture  Effect 

Further  experimental  work  done  at  ARL/PSU  119]  has  revealed 
yet  another  factor  concerning  the  influence  of  the  PV  intensity 
distribution  on  the  LSS's  measurement  accuracy.  Specifically,  it  was 
observed  that  the  PV  intensity  distribution  was  influenced  by  an 
"aperture  effect,"  which  can  be  visualized  as  follows. 

As  described  in  Section  2.1.2,  the  aperture  placed  over  the  photo¬ 
multiplier  tube  restricts  the  length  of  the  laser  beam  monitored  by 
the  measuring  subsystem  of  the  LSS;  thus,  the  aperture  defines  the 
length  of  the  PV.  Previously,  it  was  assumed  that  the  aperture 
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produced  a  sharp  definition  of  the  PV  length,  so  that  nuclei  passing 
through  the  portion  of  the  beam  monitored  by  the  system  were  always 
counted,  and  nuclei  passing  through  all  other  portions  of  the  beam 
were  never  counted.  The  "aperture  effect"  observed  at  ARL/PSU  showed 
that  these  PV  boundaries  set  by  the  aperture  were,  in  fact,  not 
sharply  defined.  Therefore,  there  is  a  dim  region  around  the  PV 
that  nuclei  can  pass  through  and  still  be  monitored  by  the  system; 
yet  because  of  the  decreased  intensity,  they  will  be  counted  by  the 
system  as  being  smaller  than  they  actually  are.  In  other  words,  there 
is  an  intensity  profile  in  the  length-wise  direction  of  the  PV  as 
well  as  in  the  width-wise  direction. 

The  impact  of  this  aperture  effect  on  the  intensity-mapping 
procedure  comes  through  step  4  of  the  analysis.  Specifically,  it 
means  that,  in  reality,  the  PV  intensity  distribution  relates  a  set 
of  normalized  PV  intensity  ranges  to  a  set  of  fractions  of  the  PV  area, 
where  this  PV  area  is  formed  by  boundaries  of  the  PV  width  and  the  PV 
length.  Thus,  the  number  of  nuclei  in  each  classification  channel  is 
given  by  the  relevant  fraction  of  the  PV  area,  not  PV  width. 

Considering  Figure  16  again,  it  can  be  seen  that  the  effect 
here  would  be  to  use  the  left-side  plot  to  become  a  two-dimensional 
intensity  distribution,  instead  of  simply  a  one-dimensional  intensity 
profile.  The  intensity  variation  in  the  length-wise  direction,  due 
to  the  aperture  effect,  would  be  similar  to  the  variation  in  the  width- 
wise  direction,  having  a  peak  intensity  near  the  center  and  dropping 
off  to  zero  at  the  extremes  of  the  PV. 


The  effect  of  such  a  two-dimensional  treatment  of  the  intensity 
variation  in  the  PV  is  predictable.  More  of  the  PV  area  will  be 
related  to  low  normalized  intensity  ranges,  and  less  of  the  PV  area 
will  be  related  to  high  normalized  intensity  ranges.  It  then  follows 
logically  from  the  intensity-mapping  procedure  that  more  nuclei  will 
be  counted  in  the  low  size  classification  channels  and  less  nuclei 
will  be  counted  in  the  high  size  categories. 

The  trend  of  this  effect  is  confirmed  by  the  results  of  the 
experiments  conducted,  as  seen  in  Table  6a.  If  the  aperture  effect 
influenced  the  counting  and  classification  of  nuclei  in  the  way 
described  above,  then  the  "Count  Rate  From  Intensity  Mapping"  (Table  6a) 
would  shift  toward  a  distribution  of  counts  that  looked  more  like  the 
"Count  Rate  From  Experiments,"  where  the  numbers  of  nuclei  counted 
are  large  in  the  lower  size  categories. 
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